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ABSTRACT 

We describe the calibration and data processing methods used to generate full-sky maps of the cosmic 
microwave background (CMB) from the first year of Wilkinson Microwave Anisotropy Probe (WMAP) 
observations. Detailed limits on residual systematic errors are assigned based largely on analyses of the 
flight data supplemented, where necessary, with results from ground tests. The data are calibrated in flight 
using the dipole modulation of the CMB due to the observatory’s motion around the Sun. This consti- 
tutes a full-beam calibration source. An iterative algorithm simultaneously fits the time-ordered data to 
obtain calibration parameters and pixelized sky map temperatures. The noise properties are determined 
by analyzing the time-ordered data with this sky signal estimate subtracted. Based on this, we apply a 
pre-whitening filter to the time-ordered data to remove a low level of 1// noise. We infer and correct 
for a small (~1%) transmission imbalance between the two sky inputs to each differential radiometer, 
and we subtract a small sidelobe correction from the 23 GHz (K band) map prior to further analysis. No 
other systematic error corrections are applied to the data . Calibration and baseline artifacts, including 
the response to environmental perturbations, are negligible. Systematic uncertainties are comparable to 
statistical uncertainties in the characterization of the beam response. Both are accounted for in the covari- 
ance matrix of the window function and are propagated to uncertainties in the final power spectrum. We 
characterize the combined upper limits to residual systematic uncertainties through the pixel covariance 
matrix. 
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instrumentation: detectors 
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1. INTRODUCTION 

The Wilkinson Microwave Anisotropy Probe (WMAP) has produced full-sky maps of the cosmic microwave 
background (CMB) of unprecedented precision and accuracy. On angular scales larger than ~ 0.°5, the dominant 
uncertainty is not the instrument noise but rather the “cosmic variance” inherent when analyzing a single realization 
(the observable universe) of a random process. In the cosmic variance limit, no further improvement can be made by 
reducing instrument noise, placing even greater importance on the minimization of non-random instrumental effects 
in the data. 

The WMAP design emphasizes control of systematic errors (Bennett et al. 2003a). The observatory was designed 
with a detailed systematic error budget in place, and a mature data analysis pipeline was written early to help inform 
many of the design decisions. Differential radiometers compare the temperature from widely-separated regions of 
the sky through back-to-back matched optics. Common-mode signals thus cancel before affecting the sky maps. 
The radiometer feed horns only illuminate a fraction of the primary mirrors, reducing the sidelobe response in the 
beam patterns. The instrument was designed to have minimal response to electrical or thermal perturbations and 
operates in an exceptionally stable environment at the second Earth-Sun Lagrange point. The observatory’s compound 
spin and precession allow rapid inter-comparison of different positions on the sky, greatly reducing the coupling of 
systematic error signals into the sky maps and effectively symmetrizing the beam response. WMAP data are calibrated 
in flight using the dipole modulation of the CMB from the observatory’s orbital motion around the Sun as a full-beam 
calibration source. We measure the beam pattern in flight using observations of the planet Jupiter. 

We characterize or limit systematic errors in the WMAP first-year data using flight data supplemented where 
necessary with results from ground tests. Systematic errors may be classified into several broad categories including 
the following: 

Calibration Errors . — The time-ordered data is simultaneously fit for the calibration and sky map. An iterative 
algorithm updates the calibration solution based on the previous iteration of the sky map solution, and vice-versa. The 
most important source of error is confusion between the dipole signal and higher-order sky signal, especially bright 
Galactic foreground emission in the low frequency WMAPbands. See §2.2.1. 

Map-making errors . — These are due to poor convergence in the sky map solution, or to errors in the determination 
of the spacecraft pointing. See §2.1.1 and §3.3. 

Beam Errors . — Instrument noise, background subtraction, and pointing errors can limit the in-flight measurement 
of the beam response from Jupiter. Although the beams are not symmetric, the observatory’s compound spin and 
precession effectively symmetrize the beam response. Uncertainties in both the beam solid angle and the window 
functions must be characterized. See Page et al. (2003a) for a complete discussion of beam mapping. We summarize 
and incorporate their results in §3.3.3. 

Sidelobe Response . — Sidelobe pickup of bright sources (e.g, the Galactic plane) introduces an additive signal 
dependent on the orientation of the beams on the sky. Barnes et al. (2003) discuss the sidelobe response of each 
radiometer and estimate the effect on the first-year sky maps. 

Baseline Errors . — Thermal or electrical perturbations can produce signals dominated by an additive term in the 
time domain. Slow drifts are removed as part of the calibration procedure, but signals near the spin period can couple 
to the sky maps with some efficiency. See §3.4. 

Striping . — Correlations in the time-ordered data from sources not fixed on the sky (e.g. 1// noise or post- 
detection filtering) introduce correlated noise in the sky maps. Application of a pre-whitening filter to the time-ordered 
data reduces this effect. See §2.3.2. 
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We have constructed a detailed model of the instrument that successfully reproduces all major aspects of the 
instrument performance. Software simulations using this model validate the map-making algorithm and allow us 
to assess the effect in the sky maps of various signals in the time-ordered data. Based on this, we apply a pre- 
whitening filter to the time-ordered data to remove a low level of 1 // noise. We infer and correct for a small (~1%) 
transmission imbalance between the two sky inputs to each differential radiometer, and we subtract a small sidelobe 
correction from the 23 GHz (K band) map prior to further analysis. No other systematic error corrections are applied 
to the data. Calibration and baseline artifacts, including the response to environmental perturbations, are negligible. 
Systematic uncertainties are comparable to statistical uncertainties in the characterization of the beam response. Both 
are accounted for in the covariance matrix of the window function and are propagated to uncertainties in the final 
power spectrum. We characterize the combined upper limits to residual systematic uncertainties through the pixel- 
pixel covariance 

This paper is organized as follows. In § 1.1 we define the terms and notation used throughout the paper. In §2 we 
discuss the iterative algorithm for making maps from time-ordered data, then generalize to the case of simultaneous 
calibration and sky map estimation. (Appendix C further generalizes to map-making with polarization data.) We also 
discuss the noise properties of the time-ordered data. In §3 we discuss combined systematic error limits from calibra- 
tion and map-making. We also present the noise properties of the sky maps in terms of their pixel-pixel covariance. 
Finally, we derive upper limits to environmental perturbations and summarize the combined systematic error budget. 
In §4 we present our conclusions. 


1.1. Notation and Overview 

Throughout this paper, we denote vectors and scalars with bold and plain lowercase letters, respectively. Matrices 
and operators are denoted with uppercase bold letters. Following Stompor et al. (2002) we denote vector and matrix 
component indices in parentheses, saving subscripts and superscripts to further identify quantities. A summary of 
the most frequently used symbols is given in Table 1. Unless otherwise stated, all temperatures are specified in 
thermodynamic units. 

WMAP measures the brightness temperature of the sky as a function of position, t(0,(f>) -4 t (/?), where p denotes 
the sky map pixel number, indexed from 0, in HEALPix nested format (G6rski et al. 1998). To make this measurement, 
WMAP scans the sky and measures the temperature difference between two points at time t. The resulting time-ordered 
differential data (TOD) is denoted d. The main goal of the map-making is to obtain the minimum variance estimate of 
the sky map, t, by inverting the raw differential data. Note that t is the true sky temperature convolved with the nominal 
instrument beam, plus instrument noise. In the process of solving for the map, we calibrate the data by estimating the 
gain and baseline from the flight data itself; characterize the full instrument beam response function from observations 
of the planets; characterize the noise spectrum of the instrument, and place limits on residual systematic errors. 

In order to produce stable data with a nearly-white noise spectrum WMAP employs 20 high-electron-mobility- 
transistor (HEMT) based differential radiometers. Each radiometer measures the brightness difference between two 
inputs, one fed by an A-side beam, the other by a B-side beam approximately 141 °apart. A detailed description of their 
design and fabrication may be found in Jarosik et al. (2003a); a summary of their in-flight performance is presented 
in Jarosik et al. (2003b). The 20 radiometers form 10 polarization-sensitive “differencing assemblies” (DA) which 
are designated based on their frequency or waveguide band: Kl, Kal, Ql, Q2, VI, V2, Wl, W2, W3, W4. The two 
radiometers in a DA are sensitive to orthogonal linear polarization modes; the radiometers are designated 1 or 2 (e.g., 
Kl 1 or K12) depending on which polarization mode is being sensed. Each of the 20 radiometers is intrinsically a 
2-channel device, with channels designated 3 and 4 in the flight telemetry, e.g., K113 or K114. [Channels 3 and 4 


-4- 


were designated left and right, respectively, in Jarosik et al. (2003a).] There are 40 such data channels in the flight 
telemetry. As discussed below, each of the 40 channels is individually calibrated, then the 4 channels from a single 
differencing assembly are combined to form differential intensity and polarization signals as follows. Let d/y be the 
calibrated differential signal from a single channel, y, of radiometer i. The differential intensity data (Stokes parameter 
/) is the average of all 4 channels 

d= ^(dl3+rii4)+-(d23+d24)- (1) 

The differential polarization data is obtained by taking the difference between the two radiometer signals 

p = -(di3 + di4)~-(d23 +^24)- (2) 

In Appendix C on polarization map making, we relate the differential polarization signal to the Stokes parameters Q 
and U. Kogut et al. (2003) discuss additional aspects of polarization mapping and analyze the first-year temperature- 
polarization correlation data based on these maps. Note that we can also form null channels from the data by taking 
channel differences, (do -d/ 4 ). As discussed in §3.2, these channel combinations provide valuable consistency tests 
for the final sky maps. 

A single channel of uncalibrated differential data may be modeled as 

c = g[M*(t+ts)+n]+b, (3) 


where each quantity is a function of time: 

c(r): raw differential data, in counts or digital units (“du”). 
g(r): instrument responsivity (here called gain), in du mK" 1 . 
b(t)‘ instrument baseline, in du. 
n(f): instrument noise, in mK. 

M * (t+ts) = At(/): differential sky signal from all sources, convolved with the instrument beam, in mK. This 
includes fixed sources, t (e.g., CMB and Galactic emission) and moving sources, t s (e.g., planets). 

In practice, the differential signal is integrated for a fixed time r and sampled at discrete times thus we may regard 
time series data as a vector with N f observations. The integration time per observation is 128.0 ms, 128.0 ms, 102.4 
ms, 76.8 ms, and 51.2 ms for bands K through W, respectively. 

The differential temperature at time t is the convolution of a time-dependent mapping function, M, with the sky 
signal at time t , t+t s (r) 

At(f) = J dCl nl M(n(f ) , n') (t(u') + ts(n' , t)). (4) 

Here t = tc + t g is the fixed sky signal consisting of CMB anisotropy, L, and Galactic foreground signals, t g , while t* 
represents all time dependent sources, especially the Sun, Earth, and Moon, which are potentially visible in the far 
sidelobes of the instrument. The operator M can be represented as an N T x N p matrix where each row is the normalized, 
full-sky beam response in sky-fixed coordinates as given by the scan pattern. Several features of the mapping function 
that pertain to the data processing are discussed in Appendix A. The main beam response is mapped in flight using 
observations of Jupiter as a far-field point source (Page et al. 2003a). An important aspect of the WMAP optical design 
(Page et al. 2003b) was to limit sidelobe pickup to negligible levels and to have the effective beam response in the final 
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sky maps be approximately circularly symmetric. We discuss each of these topics in more detail in separate papers, 
(Barnes et al. 2003; Page et ai. 2003a), while this paper summarizes the main results in terms of systematic error limits 
in the sky maps. 


The instrument gain, baseline, and noise are determined from the flight data itself. This is an iterative process that 
we discuss in detail in §2.2.1. Here we provide a brief overview of our terminology in order to frame the following 
discussion of systematic errors. Let the gain and baseline measured in flight be g and b, respectively. The calibrated 
differential signal is then 


a= (c-« = ! i . t+ i ii+ b-b 

g g g g 

With calibrated data available, the sky map is obtained by evaluating the linear equation 


(5) 


t = Wd, 


( 6 ) 


where W is a linear operator defined in §2.1. The properties of W are determined by the scan strategy of the observa- 
tory and the noise properties of the time-ordered data, d. 


The WMAP scan pattern is an integral part of the mission design (Bennett et al. 2003a). It consists of a compound 
spin and precession centered about the Sun- WMAP line, with parameters as given in Table 2. There are several 
aspects of this scan strategy that are important for high quality data: scans of a given pixel cross at many angles so 
that the effective beam response is symmetric; a given pixel is observed on many different time scales from minutes 
to months; the angular velocity of a given line of sight is nearly constant on the sky; the instrument observes a large 
fraction (>30%) of the sky each day; and the time-average of the differential data is approximately zero over an hourly 
calibration period, allowing for robust initial baseline estimation. 


2. THE MAP-MAKING PIPELINE 

A graphical overview of the WMAP data processing and analysis pipeline is shown in Figure 1 of Bennett et al. 
(2003b). The heart of the pipeline is a set of programs that bring science and housekeeping data from the Science and 
Mission Operations Center (SMOC) through to a set of calibrated full sky maps for each of the 10 WMAP differencing 
assemblies. Numerous additional programs are used to generate ancillary data products such as beam response maps, 
calibration files and analysis products. In this section we provide a high-level overview of each program in the map- 
making pipeline and provide references to later sections of this paper, or to companion papers, for further details. 

Raw telemetry data from the satellite is transferred approximately once per day through NASA’s Deep Space 
Network (DSN), to the SMOC, located at the Goddard Space Flight Center. The SMOC monitors the basic health and 
safety of the Observatory, sends all commands, and requests re-transmissions of data that were missing from a previous 
transmission. The data are then “level-0” processed into a set of time-ordered, daily files which contain science data, 
instrument housekeeping data, spacecraft data (including attitude and ephemeris data), and event message files. These 
files are then transferred to the Science Team’s processing facility, also at the Goddard Space Flight Center. 

Every time a new full day of data arrives, a series of automated procedures perform the following tasks: 1) Gen- 
erate a standard set of daily plots that are archived and visually inspected; 2) Generate a reduced “trending archive”, 
which consists of subsets of the data sampled once every 10 minutes. In the case of the science data, we record the 
mean and rms of each channel in a 10 minute interval, whereas the housekeeping data are sub-sampled, once per 10 
minutes; 3) Perform a series of data quality checks that search the data for violation of pre-set range limits or excessive 
time-gradients in the telemetry signals. Limit violations are logged and notification is sent to a member of the science 
team via e-mail. (Initial limit tests are performed at the SMOC as well.) 
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At selected time intervals, the level-0 telemetry files are collated by a pre-preprocessor into a master archive of 
raw (uncalibrated) data. The major functions of the pre-processor are to: 1) collate the science and housekeeping data 
into single daily files; 2) flag data that is suspected or known to be unusable; 3) interpolate the attitude and ephemeris 
data to times that are commensurate with the science data time stamps; and 4) apply a coarse flag to data that is within 
7° of one of the outer planets (Mars through Neptune) so that it may be rejected from the initial sky maps, but identified 
for beam mapping (Page et al. 2003a). (Only Jupiter data is used for making the final beam maps.) 

Initial sky maps and calibration data are generated from the raw archive using the iterative map-making algorithm 
first described in Wright et al. (1996), and further described in §2.1. As discussed in §2.2, the initial calibration is 
determined by fitting the raw time-ordered data to the known signal produced by the CMB dipole. Because the sky 
signal contains significant higher-order power (/ > 1), the calibration solution must be iteratively improved in concert 
with the initial sky map iterations. The convergence of this simultaneous fit has been demonstrated with end-to-end 
simulations, which are also described in §2.2. 

The calibration solution converges more rapidly than the sky map does, so we freeze the initial calibration solution 
after ~ 10 iterations before proceeding to convergence with the sky maps. With the initial dipole-based calibration data 
in place, a re-processor generates a refined gain and baseline solution and applies this to the data. The program also 
updates the data quality flags, as necessary, then writes the calibrated data to a new final time-ordered archive. The 
refinements to the dipole-based gain solution are discussed in Jarosik et al. (2003b) and in §2.2.1. The initial baseline 
solution is refined with a pre- whitening filter (Wright 1996) which is presented in detail in §2.3.2. 

The final sky maps are computed using the final calibrated time-ordered data as input. The first-year sky maps 
required 20 post-calibration iterations to be sufficiently converged. The map-making algorithm is fundamentally the 
same as is used in the initial estimates, but we add some refinements for this final stage of processing. These include: 1 ) 
correcting for a small ( £ 1%) loss imbalance between the A and B-side sky beams. Jarosik et al. (2003b) demonstrate 
that this effect is nearly orthogonal to the gain solution, so its inclusion after the calibration processing does not 
invalidate the gain solution. 2) We weight individual time-ordered observations by their proper statistical weight to 
account for the small change in instrument noise (< 1 %) over the course of a year due to the 0.9 K seasonal temperature 
variation experienced by the instrument cold stage. 3) We compute the planet-boresight angle for each observation to 
minimize unnecessarily conservative data loss. The criterion used for the first-year maps is a cut of radius 1 .°5. 

We also generate polarization maps using a generalization of the temperature algorithm, the details of which 
are presented in Appendix C. For the first-year data processing, we have generated maps of the Stokes parameters 
Q and U, but we have not yet fully characterized all of the potential systematic errors in these maps. However, 
the temperature-polarization correlation data are much less prone to systematic errors than the polarization auto- 
correlation data. Kogut et al. (2003) have analyzed the temperature-polarization data and, supported by systematic 
error limits from Barnes et al. (2003), they find a significant correlation, including the signature of a relatively early 
epoch of cosmic reionization. 

In the remainder of this paper, we present a more detailed description of each of the map-making and calibration 
procedures, including an assessment of their performance with the first-year WMAP data. We then derive detailed 
systematic error limits applicable to the time-ordered data, to the sky maps, and to the angular power spectrum. 


2.1. Map-Making with Pre-Calibrated Data 


While the process of generating the final sky maps from calibrated data comes last in the map-making pipeline, 
we discuss the algorithm first because the algebra of map-making is central to the entire data processing scheme, and 
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it helps to guide the systematic error analysis. 

We consider the problem of estimating a sky map, t, from calibrated, differential time-ordered data, d, which is a 
linear function of the sky map 

d = Mt+n, (7) 

where M is the mapping function of the experiment, which has N p columns and N r rows, N p is the number of sky map 
pixels and N t is the total number of time-ordered observations. In its simplest form, each row (observation) of the scan 
matrix contains a +1 in the column (pixel) seen by the A-side beam, and a -1 in the column (pixel) seen by the B-side 
beam. This matrix can be generalized to include the effects of beam convolution, but for WMAP these refinements 
are small and are being deferred to future processing. The effects of a differential loss imbalance between the A and 
B-side beams is readily accounted for by using values different from ± 1 in M. An analysis of this effect in the WMAP 
radiometers is presented in Jarosik et al. (2003b). The details of how we account for it in the pipeline are given in 
Appendices A and C. 

The noise n is assumed to have zero mean and covariance N, 

(n) = 0 (8) 

<nn r ) = N. (9) 

We defer a detailed discussion of the WMAP noise properties to §2.3.2, but for most radiometers it is reasonable to 
approximate the noise covariance as diagonal, N ~ <Tq I, where ao is the rms noise per observation and I is the identity 
matrix, though this assumption is not required for the algorithm described below to converge. 

The least-squares, or maximum-likelihood estimate of the sky map, t, results from solving the normal equations 

t = (M r N _1 M) _1 -(M r N -1 d). (10) 

More generally, we obtain an unbiased estimate of the sky map by choosing any symmetric matrix S in place of N _l 
(Tegmark 1997). To see this, substitute Mt+n in place of d in the above equation to get 

t = (M r SMr‘ •(M r S[Mt+n]) = t+(M r SMr 1 -(M r Sn). (11) 

Thus t reduces to t plus a noise term that is independent of t and has zero mean over an ensemble average. For the 
first-year processing we take S = I. To simplify notation, we define a matrix W = (M r M) _1 • M r in which case the 
map solution is t = Wd. 

The pixel-pixel noise covariance in the sky map solution is 

£ = <(t— t)(t-t) r ) = ((Wn)(Wn) r ) = WNW r 

= (M r M)-' • (M r NM) • (M r M) _ 1 . (12) 

In the limit that N is diagonal and the rms noise per observation, ao, is constant, £ reduces to ao(M r M) -1 . For 
the WMAP scan pattern, the matrix M r M is diagonally dominant with diagonal elements n 0 b S (p)> the number of 
observations of pixel p by either the A or B-side beam. Thus, to a very good approximation, the pixel-pixel covariance 
matrix is diagonal 

Values for (tq are given by Bennett et al. (2003b). 

The leading order off-diagonal terms occur at the beam separation angle (0beam ~ 141°), and are of order 0.3% of 
the diagonal elements. If the time-ordered noise N is not diagonal, then maps produced with the above algorithm will 
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have correlated noise (stripes) along the scan paths defined by M. This is a small, but not negligible, effect for some 
of the WMAP radiometers, and is entirely negligible for others. The noise properties of the time-ordered data and sky 
maps are further discussed in §2.3.2 and §3.2, respectively. 


2.7.7. Iterative Map Making 


The evaluation of the sky map solution Wd requires the inversion of the N p x N p matrix D = M r M. We use the 
iterative approach introduced by Wright et al. (1996) to evaluate this expression. Briefly, suppose we have an initial 
guess for the sky map, to, which differs from the true sky map, t, by Jto = to - 1 . Then D to = D (t + 6to) can be recast as 

D<5to = Dto-M r d, (14) 


where we have used the fact that Dt= M r d. As noted above, D is diagonally dominant, so a good approximate inverse 
for D is 


6 '(PhPj) 


nobs( Pi) ‘ J 


(15) 


This leads to the approximate solution for the residual, Sio ~ D ^Dto-M^d], and suggests the following iterative 
solution 



06) 

t„-b-‘[Dt n -M r d] 

(17) 

(D -1 M r )d + (I-D _1 D)t„. 

(18) 


The interpretation of equation (1 8) is that for each pixel the new sky map temperature is the average of all differential 
observations of that pixel (accounting for the sign of the observing beam) corrected by an estimate of the signal 
in the paired beam, based on the previous sky map iteration. The expression in equation (18) can be efficiently 
evaluated because the sums can be accumulated by reading through the time-ordered data from disk, each iteration, 
and accumulating data into arrays of length N p 


n 0 bs * tn+\(PA) ~ ► n 0 bs ■ t re+ i (pA) + Wi [dfeHt^p#)] n 0 bs (Pa) n obs(/?4) + w/ 

n 0 bs * Wp B ) n 0 bs 't n +\(pB)-Wi [&(ti)-t n (pA)] n obs (pb) n obs (p B )+Wi, (19) 

where wi = 1 in the initial sky map processing and is proportional to a noise weight (equation C31) in the final sky 

map processing. Note that it is never necessary to store or invert an N p x N p matrix. 

We have tested this algorithm extensively with flight-like simulations. In this section we present results for an 
“ideal” noiseless instrument with circular beams and perfect calibration to isolate the performance of the map-making 
algorithm from other effects. More realistic data models are introduced to the simulation in subsequent sections. 
Figure 1 shows a sample residual map, t our -ti n , generated from a one-year simulation of Q2 data. The input sky 
map included realistic CMB signal with a peak-to-peak amplitude of ~ ±420 pK, and a Galactic signal with a peak 
brightness of ~50 mK. The output sky map is recovered with an rms error of < 0.2 /zK, after 50 iterations. The 
dominant structure in the residual map is a mode aligned with the ecliptic plane. The power in this mode is concentrated 
in spherical harmonic mode / = 4, due to a combination of the WMAP scan strategy and the beam separation angle. 
This is the mode on the sky that is least well measured by WMAP (except for the monopole!) and is thus the slowest 
to converge, though additional iterations would reduce its amplitude even further. The final first-year flight sky maps 
were effectively run for 80 iterations. Since the rms error associated with this term is very small, and since we build 
up a more realistic data model in subsequent simulations, we do not further quantify this contribution to the final 
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systematic error budget. Rather, we subsume it into an overall map-making and calibration error budget that includes 
this and other effects together. 

This iterative approach to map-making is readily generalized to polarization maps as well - the formalism is 
presented in Appendix C. We have tested that algorithm with the same simulations used above to test the temperature 
algorithm and find that the polarization maps converge even faster than the temperature maps. After 10 iterations, the 
map-making artifacts in a residual polarization map are <0.05 p K peak-peak. 


2.2. Combined Calibration and Map Making 

The processing algorithm described above assumes that the data have already been calibrated. In practice, we use 
the above algorithm in the second stage of map-making, after an initial stage in which we simultaneously solve for the 
radiometer calibration and the sky map. In the initial stage of map-making, we employ the same iterative algorithm 
to solve for the map, but rather than processing straight through the time-ordered data on each iteration, we process 
the data one hour at a time, pausing to solve for the calibration in each radiometer channel, before accumulating the 
calibrated data. The calibration solution then iteratively improves as the sky model improves. The following sub- 
sections lay out the procedure in detail, and present results for the flight data with an assessment of its precision and 
accuracy based on flight-like simulations. 


2.2. 7. Instrument Calibration from the Dipole Modulation 

For a sufficiently short period of time the instrument gain and baseline can be approximated as constant, c ~ 
g*(At+n) + fc*, where g k and b k are the gain and baseline during the k tb 1-hr calibration period. Since the sky signal 
At is dominated by the CMB dipole measured by COBE> Atj, a single channel of raw data can be modeled as 

CmiSkibfc) = £jt(At<j + At v ) + 7ty, (20) 

where At v is the additional dipole moment induced by the motion of WMAP relative to the solar system barycenter 
(the rest frame of the COBE dipole). 

We fit for the gain and baseline in each calibration period, k , by minimizing 

■^2 . [cfr)-C m (f/|g*A)] 2 (21) 

/€* 

where i is a time-ordered datum index. We omit data that are flagged as unusable, and data when either the A or B-side 
beam points within a Galactic pixel mask. The mask used for this latter application is the Kp8 mask defined in Bennett 
et al. (2003c), without edge smoothing. This mask is used throughout the map-making pipeline. The fit is performed 
for each of the 40 WMAP channels independently. To minimize the covariance between the recovered gain and the 
baseline, it is necessary to have a scan strategy such that the time average of the sky signal. At, is approximately 
zero in one calibration period. The combined spin and precession of WMAP is designed to produce time-ordered data 
that satisfies this requirement. For example, in K band, which has the largest sky signal, a 1-hr running mean of the 
differential sky signal has an rms fluctuation of 14 pK, compared to a dipole signal of greater than 3 mK. After each 
hour of data is processed for the calibration solution, the data are accumulated as per equation (19) to develop the sky 
map solution. 
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The largest source of error in the calibration fit is due to un-modeled sky signal from the CMB anisotropy and 
Galactic foreground emission, At a = At- At d . This projects onto the dipole signal and, as shown below, causes errors 
in the gain solution as large as 5-10% in K band, where the Galactic signal is strongest. The calibration fit may be 
iteratively improved by subtracting an estimate of the anisotropy from the raw data prior to fitting . In particular, let g[ 
be the gain inferred for calibration period k from the previous iteration of the calibration fit. Then minimize 


* 2= E 

i£k 


_2 5 


( 22 ) 


where 

c' = c-g'A< (23) 

and At* is the differential sky signal (less the dipole component) computed from the previous sky map iteration. This 
process is repeated until the calibration solution is sufficiently converged. 


Note that the absolute calibration is tied to the time-dependent portion of the dipole signal, At v ; we use the fixed 
dipole as a short-term transfer standard only. In particular, when we update the sky model and apply the anisotropy 
correction in equation (23), any error in the fixed dipole moment, At d , is assumed to be anisotropy, and is applied as a 
correction in the same w ay. For a data set of at least one year in length (one full cycle of At v ), the error in the absolute 
calibration will be essentially orthogonal to any error in the fixed dipole At d . 


2. 2. 2. Performance of the Dipole-Based Gain Solution 

As an illustration of the systematic gain error induced by higher-order (/ > 1) anisotropy, Figure 2 shows an 
example of the gain solution convergence from a one-year low-noise simulation. This simulation implements the 
simultaneous calibration and sky map estimation discussed above and was run for 30 iterations. The example shown 
is for one channel of K band data (the worst case) which exhibits ~7% errors after one iteration, corresponding to a 
sky model that has only a dipole component. After 30 iterations, the residual errors are <0.1 % over the entire year. 
Similar, or better, performance is achieved for all other WMAP channels. 

In processing the final first-year maps, the combined calibration and map-making code was run for 10 iterations. 
However the initial sky model was based on an earlier “pathfinder” run of the pipeline that ran for a total of 30 
iterations of combined calibration and map-making plus an additional 20 iterations of sky map convergence. Thus 
we conservatively estimate that the combined absolute and relative calibration errors due to incomplete calibration 
convergence to be <0.1%. We defer a discussion of the final calibration uncertainty to §3.1. 

Figure 3 shows a sample of the converged gain solution from equation (22) for two WMAP channels, K1 13 and 
VI 13. Note that the radiometer gains are typically drifting by a few percent over the course of the first year. As we 
show below, the dipole-based fits easily track drifts at this level. The noise in the gain solution is typically a few 
percent per hourly calibration period, though, as is readily seen in the figure, the noise level changes with time of 
year as the scan pattern sweeps around the CMB dipole. The VI 13 gain exhibits an additional modulation that is 
clearly correlated with the physical temperature of the instrument. However, the time scale of the temperature change 
is slow enough that the corresponding gain changes are well tracked by the dipole fits. Quantitative limits on thermally 
induced gain and baseline errors are discussed in Jarosik et al. (2003b) and in §3.4.1. A summary of the gain statistics 
from the flight data is given in Table 3. 


I 
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2.23, The Initial Baseline Solution 

The bottom two panels of Figure 3 show the converged baseline solution resulting from the fit in equation (22) 
for one year of K1 13 and VI 13 flight data. The fits have had a mean subtracted, and have been divided by the 
gain to convert to temperature. These plots, which are representative of all 40 channels, show that the offsets of the 
radiometers are typically stable to ±5 mK over the course of the first year. Simulations demonstrate that the hourly 
baseline solution is unbiased. However, it is also clearly noisier than optimal, consistent with the flight measurements 
of the noise power spectral density (Jarosik et al. 2003b). In §2.3.2 we describe an improved baseline model that is 
based on the application of a pre-whitening filter tailored to the measured noise spectrum of each channel. 


2.3. Improving the Calibration Model 

The sky maps obtained with the hourly calibration are reasonable; however the noise in the calibration solution, 
particularly in the baseline, is significantly higher than optimal, and the use a piecewise continuous calibration in the 
final maps would introduce striping in its own right. In the following subsections, we present the steps undertaken to 
filter the gain and baseline solutions that enter into the final sky maps. Prior to generating the final maps, this refined 
calibration is applied to the data and written to disk as the final first year calibrated time-ordered archive. 


23.1. The Gain Model 


Jarosik et al. (2003b) present a physical model for the gain that is based on the RF bias, or “total power” measured 
in each channel, and on the physical temperature of the instrument cold stage, which is monitored with high resolu- 
tion platinum resistance thermistors (PRTs). Each of these quantities is recorded once every 23 s in the engineering 
telemetry with a relative noise that is substantially lower than the noise in the dipole-based gain solution. Thus, if the 
model fits the dipole-based data satisfactorily, it offers a means for measuring the gain with more precision, and on 
time scales shorter than the spin period. The model for the gain, g(r), has the form 


g (O = go 


m-v o 
wo-iv 


(24) 


where V is the measured RF detector bias, 7}^ is the measured temperature of the FPA, and go, Vo, and To are fit 
constants. See Jarosik et al. (2003b) for more detail. 


Figure 8 in Jarosik et al. (2003b) shows the performance of the gain model when fit to the dipole-based gain 
solution (see also Figure 7 in this paper) . In § 3 . 1 we evaluate the overall performance of the gain model and the hourly 
gain solution in the context of an end-end simulation designed to place limits on combined calibration and map-making 
errors. We will conclude that the model provides an excellent description of the radiometer gain, and here we adopt it 
as the final gain solution for further processing. The gain model fits into the data processing sequence as follows. After 
we iterate the simultaneous calibration and sky map solution long enough for the calibration to converge (10 iterations 
when starting with a good sky model), we freeze the dipole-based calibration and fit the gain model parameters in 
equation (24). All subsequent data products are produced with data calibrated using this gain solution, including the 
time-ordered archive, the final sky maps, and the Jupiter beam maps. 
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2.3.2. Baseline Filtering 

The baseline that results from the initial calibration solution is not optimal. This is due to the fact that the initial 
baseline is sampled once per hour (0.28 mHz), while Jarosik et al. (2003b) show that the power spectral density of the 
noise has a 1// knee frequency of a few mHz, typically. If the initial baseline estimate were used in the final sky maps, 
it would generate weak stripes of correlated noise along the scan paths, as per equation (12). Even so, it is important 
to note that 1 // effects are small relative to the white noise. In the worst WMAP radiometer, W41, the amplitude 
of the noise covariance N at small lag is ~2% of the white noise variance. Thus we treat 1 // noise iteratively in the 
data processing by first ignoring it to obtain an estimate of the gain, baseline, and sky solution. Then we subtract the 
estimated sky signal from the time-ordered data, apply a pre-whitening baseline filter to the residual noise, add the 
sky signal back in, and write the data to a final calibrated, time-ordered data archive. The approach of first subtracting 
an estimated sky signal is designed to avoid biasing the gain solution and/or removing low-order power from the sky 
maps. The noise properties of maps constructed in this way must account for the filtering process. We discuss the 
map-making algebra appropriate to our filter implementation in Appendix B. 

The steps we follow to define and apply the filter are as follows. We remove an estimate of the sky signal, in du, 
from the raw differential data using 

c'(0 = c(f)-g(f)At'(r)-fc* ? (25) 

where b k is the hourly baseline point appropriate to the current time, g is the final gain solution from equation (24) and 
At' is the differential sky signal computed from the initial sky map. We then evaluate the auto-correlation function of 
c'(f) to a lag of 10 4 sec. Results for representative radiometers are shown in Figure 4. This range of lags is sufficient 
to account for both the long-range correlations due to 1// noise and the correlation at a lag of 1 observation due to the 
low-pass post-demodulation filter in the Analog Electronics Unit (Bennett et al. 2003a). The baseline filters are then 
defined as follows (Wright 1996) 

1. Fit the auto-correlation function, C(Ar), to the model defined below. 

2. Fourier transform the model correlation function to generate the power spectral density P(f). 

3. Compute the Fourier space filter w(f) = 1 />/P(f) and set w(0) = 0 to produce a zero mean output signal. 

4. Fourier transform the filter w(f) to generate the time domain filter w{i), normalized to 1 at / = 0. 

By inspection, the auto-correlation functions are well modeled by a log-linear function 

{ Ci A t-r 

A-£log(A*/l s) r < At < 10*/* s (26) 

0 Ar > 10 A t B Sj 

where r is the integration time for a single observation, Ci is the correlation at lag r, measured from the data, and 
A and B are fit parameters. Note that A gives the typical fractional covariance at small lag, while the suppression of 
correlations at large lag (~2000 s) is dictated by the subtraction of the hourly baseline as a pre-filter. The best-fit 
parameters are given in Table 4, and fits for selected radiometers are shown in Figure 4. 

The derived pre-whitening filters, w(f), are plotted as a function of frequency in Figure 4. One point of particular 
interest is the filter response at the spin frequency, 7.7 mHz. As shown in Table 4, the best channels have a transmission 
of > 95%, while the worst case, W41, is just above 35%. These values indicate the amount by which the dipole 
(calibration) signal would be suppressed if the filter were applied prior to calibration and sky signal subtraction. The 
convolution of c l (t) with w(t) is performed in Fourier space using the Numerical Recipes routine convlv (Press et al. 
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1992). The number of data points convolved at any one time is chosen to be the smallest power of 2 such that the data 
span a full day with sufficient padding beyond the day to guarantee that wrap-around effects are negligible. This is 2 20 
for K-Q bands, and 2 21 for V,W bands, which gives a minimum of 2.9 hours of padding on each end of a day. Sample 
auto-correlation functions obtained from the filtered data are shown in Figure 4. The filtering is clearly effective at 
removing low frequency noise in the time-ordered data. Another example of filtered data is seen in Figure 5, which 
shows 1 day of W42 data, one of the worst radiometers for 1 // noise, before and after filtering. These data are 
smoothed with a 46 s window to show structure in the unfiltered data since plots of unsmoothed data before and after 
filtering are virtually indistinguishable. 

The above results are encouraging but not definitive, because the process of sky signal subtraction and re-addition 
could introduce correlated artifacts that these tests would miss. The ultimate test of a filter is its ability to “clean” the 
pixel-to-pixel covariance matrix of the final sky maps and the noise covariance of the angular power spectrum, without 
altering the underlying sky signal. The sky map noise properties are discussed in §3.2, while the noise properties of 
the power spectrum are quantified in Hinshaw et al. (2003b). 


2.33 . Baseline Jumps 

Limon et al. (2003) identify 21 instances of sudden baseline jumps, or “glitches”, during WMAFs first year 
of operation. These events have been identified as small shifts in the properties of several microwave components 
resulting from sudden releases of internal mechanical stress, presumably from thermal changes. These events last for 
less than 1 s, and cause no discernible change in the radiometer gain or noise properties. 

Care must be taken in the application of the baseline filter in the vicinity of these steps to avoid ringing in the 
filtered data. Each event is initially flagged by the pre-processor for at least 1.2 hours on either side of the event. 
Since the initial hourly baseline is derived entirely in the time domain, the ±1.2 hr flagged interval ensures that this 
baseline estimator only “knows” about data on one side of the jump or the other. Prior to convolving the raw data 
with the pre-whitening filter, we subtract the initial hourly baseline from the data as a pre-conditioner. Thus all data 
that is input to the convolution routine has approximately zero mean. On output, the re-processor expands the flagged 
interval by 0.5 hr on either side of the event to ensure that no edge effects propagate into the usable data. In the first 
year of operation, a total of 0. 13% of the data was lost to these steps. See Table 2 in Bennett et al. (2003b). 

The threshold amplitude for jump detection by visual inspection is ~0.05 du, which corresponds to a jump of 
~150 pK in the calibrated output of radiometer W 12, the worst offender. To assess the effect of undetected baseline 
steps in the data, we have generated a test data set in which we take 24 hours of flight W12 data and insert a step of 
0.05 du in each channel. We then run the data through the pre-whitening filter to see the magnitude of the transient 
response. The result is a transient baseline error with a peak magnitude of 80 pK, which lasts for less than one 2- 
minute spin period. The total time the baseline error exceeds 10 is 22 minutes, or approximately 1 1 spin periods. 
We pessimistically assume that there could be as many as 40 steps at or just below the threshold of detection, and that 
half of these are in W 1 2. If we assume these occur at random times, and note that WMAP observes ~30% of the sky in 
any given hour, then any given sky pixel is likely to “see” approximately 1 1 x 20 x 0.3 66 data points with baseline 

errors greater than 10 zzK. Since the sign of a given step is random, and since W1 1 data is combined with W12 in the 
sky maps, we estimate the residual systematic error in a given pixel of the W1 map is less than 10 pK /2/\/20 ~ 1 
pK. We emphasize that no jumps have been observed in any other W band radiometer, thus DA -DA consistency is an 
excellent test of whether any statistic is sensitive to baseline errors of this nature. We have found no evidence that the 
W1 map is “out of family” (§3.2). 
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2.4. Final Sky Map Processing 

Once the calibrated, time-ordered archive has been written, final sky map processing commences based on the 
algorithm presented in §2.1. At this stage in the sky map processing, we add a few features to the algorithm that, 
for simplicity, are not present in the combined calibration and map-making code. These include: 1) Weighting each 
datum by a true weight, 1 /of based on an estimate derived from the physical temperature of the instrument cold stage. 
This introduces ~1% variations in the data weights over the year, since the instrument noise is a weak function of 
temperature, and the temperature varies by ~1% over the year. 2) Accounting for loss imbalance, as discussed in 
Appendix A. In effect, we model each differential observation as At = (1 +x im )t(pA)-(l-Xi m )t(pB), where x,* m is the 
small loss imbalance parameter given by Jarosik et al. (2003b). 3) Computing the planet avoidance flag at run time to 
reduce the amount of data lost. In the final sky maps, a total of 0. 1 1 % of the data was lost to planet avoidance. See 
Table 2 in Bennett et al. (2003b). 

The final stage of sky map processing, based on the filtered data, was run for a total of 20 iterations. Convergence 
was determined by measuring the rms difference between pairs of iterations for a given differencing assembly. For 
example, the difference between the 10th and 20th iteration of the W2 sky map is 0.08 /xK rms. We estimate that 
artifacts due to lack of solution convergence are <0.1 /xK rms with all of the power being in the low multipoles, 
/ < 10. We present a final combined estimate of sky map artifacts due to calibration and map-making errors in §3.1. 
This estimate includes the convergence limits given above. 

To assess the effect of the improved calibration model on the final sky maps, we form differences between the 
final post-filtered maps and the last iteration of the intial, pre-filtered maps. The results for DAs W3 and W4 are shown 
in Figure 6. The top panels in this Figure show the difference maps from a one-year simulation (§3.1) that included a 
realistic radiometer noise and gain model. The bottom two panels show differences from the flight data. Because these 
maps are largely based on the same data, most of the white noise drops out of these differences. The remaining “blobs” 
of white noise result from the change in the planet cut and can be ignored. The striking feature is the striping present 
in the W4 difference, but virtually absent in the W3 difference. As we show in Appendix B, this is the structure that 
has been removed from the data by the pre-whitening filter, an interpretation that is substantiated by the analysis of 
the simulation. The fact that the W3 difference is very small is an indication that the level of striping in the unfiltered 
W3 data was very small to start with. We estimate the level of residual striping in the final maps in §3.2. Images of 
the final maps at each frequency are presented by Bennett et al. (2003b). 


3. SYSTEMATIC ERROR ANALYSIS 

As discussed in the introduction, systematic errors may be classified by the nature of their source. In this section 
we place limits on the level of systematic errors in the final sky maps, using that classification to guide the analysis. 
In §3.1 we place limits on combined calibration and map-making artifacts, based largely on a detailed simulation of 
the first year of WMAP operation. In §3.2 we present null tests based on difference maps formed from a variety 
of data combinations, each of which should yield no sky signal. We use these maps to measure or place limits on 
correlated pixel noise (striping) in the final first-year maps. In §3.3 we discuss systematic errors relating to pointing 
and beam mapping errors. We conclude by placing stringent limits on residual errors due to environmental (thermal 
and electrical) and other miscellaneous sources. 
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3.1. Calibration and Map-Making Errors 

To assess the combined errors from calibration and map-making artifacts, we have generated a high fidelity 
simulation that includes all of the effects we believe are important for calibration and map-making. In particular this 
simulation includes: 1) A sky model that closely mimics the statistical properties of the observed sky; 2) A realistic 
noise model for every channel, including 1 / / noise [see Jarosik et al. (2003b) for a tabulation of 1 // knee frequencies]; 
3) A model for the thermal drift of the gain, baseline and offset of each radiometer, based on measured susceptibility 
coefficients, and driven by the actual temperature profile measured in flight. This simulation generates the sky signal 
using a circular beam approximation. The effects of elliptical beams are treated in a separate, noiseless simulation in 
§3.3.4. We write simulated science data to files that mimic the raw telemetry, then process the data using the same 
pipeline as was used to process the flight data. 

The top panel of Figure 7 shows the converged gain solution from the simulation for channel VI 13; the bottom 
panel shows the corresponding result from the first year of flight data. In both panels, the “noisy” black traces are the 
hourly gain data, binned in 24-hr samples to reduce the noise, and the green traces are the best-fit gain model (§2.3.1). 
For the simulation, the input gain used to generate the data is shown in grey. The absolute gain is recovered in the 
simulation to better than 0.1% for all 40 channels. 

The dipole signal seen by an observer moving with speed v relative to the rest frame of the CMB is To v/c, where 
To is the absolute temperature of the CMB, and c is the speed of light. Thus, additional sources of error that could 
affect the absolute calibration of the WMAP data include errors in the determination of WMAP s velocity with respect 
to the solar system barycenter (the point of reference for the COBE dipole) and errors in the absolute temperature of 
the CMB. The velocity of WMAP is routinely measured with respect to geocentric inertial coordinates (GCI) with an 
accuracy of <1 cm s" 1 . The velocity of the Earth is determined from the JPL ephemeris with similar accuracy. The 
combined uncertainty from velocity errors is 0.1 nK. The uncertainty in the absolute temperature of the CMB is 0.1% 
(Mather et al. 1999). Combining these uncertainties with the results of the simulation, we conservatively estimate an 
absolute calibration error of 0.5% for the first-year WMAP data. 

Errors in relative calibration can produce structure in the sky maps, beyond an overall normalization factor. The 
largest relative discrepancy between the dipole gain solution and the gain model in the flight data is ~0.4% in K 
band, and ~0.2% in the other bands. Similar deviations are seen in the simulation, thus we use the simulation as our 
primary tool for placing systematic error limits due to relative calibration and map-making errors. We have generated 
residual maps from the simulated data by subtracting the known input sky signal from the maps produced by the 
pipeline. These residual maps exhibit no visible structure aside from the pixel noise. In order to assess the errors due 
to map-making artifacts, we compute the angular power spectrum, C/, of the residual maps and search for features 
in the spectra beyond a simple flat, white noise spectrum. The residual spectra for all 10 differencing assemblies are 
shown in Figure 8 and summarized in Table 5. In general, the spectra are consistent with white noise over a wide 
range of multipole moments, but clearly show the most variation at low /. Because of this, we specifically highlight 
these modes in Table 5, where we giveC 2 , (C/) 3 _ 10 , and (Q) u _ m for each of theDA’s. For combined systematic error 
limits due to calibration and map-making, we assign twice the excess variance in each / range relative to the white 
noise plateau, a sys = 2 |(Q) b and~^) 700 -iooo|- These values are also quoted in Table 5. For comparison, the average 
power in the CMB in each band is C 2 ~ 130 /iK 2 , (C/) 3 _ I0 ~ 150 fi K 2 , and (C/) 11 _ ]00 ~ 6 fiK 2 . 

Because the simulation includes realistic models of 1 // noise and long-term thermal effects, these limits also 
implicitly limit artifacts at low l due to these effects. As we demonstrate in subsequent section, we feel this simulation 
captures all of the important radiometric characteristics of the instrument. Potential artifacts due to optical effects, 
especially pickup through the far sidelobes, are treated in Barnes et al. (2003), and are summarized in §3.3.5. 
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3.2. Difference Maps and Noise Properties 

Difference maps are combinations of the data that, ideally, should contain no sky signal. They provide insight 
to potential systematic errors and can be used to characterize the noise properties of the sky maps. The first set of 
difference maps we generate are between DA pairs with the same frequency and beam response, namely ^(Ql-Q2), 
|(V1— V2), and £(W12-W34), where W12 = i(Wl+W2), and W34 = ±(W3+W4). Images of these difference maps 
are shown in Figure 9, along with low resolution versions of the sum (signal) maps to give a sense of the signal 
strength in each map. Aside from the pattern of the noise, which follows the sky coverage [see Figure 3 of Bennett 
et al. (2003b)], the only visible structure in these difference maps is in the Galactic plane, especially in V band. This 
is understood to be a result of a small difference in the effective center frequency of the VI and V2 differencing 
assemblies (Jarosik et al. 2003a). In particular, the VI map has an effective frequency approximately 1 GHz lower 
than V2. Since the spectrum of the Galaxy at V band follows T A (v) ~ v~ 2 (Bennett et al. 2003c), we expect the Galactic 
signal to be ~3% brighter in VI than V2, which is consistent with the residual signal seen in the difference map. [A 
complete tabulation of effective center frequencies, radiometer by radiometer, is given by Jarosik et al. (2003a) for 
diffuse sources, and by Page et al. (2003a) for point sources.] Note that because the data are calibrated using the CMB 
dipole, there should be no residual CMB signal in such a difference map. A more sensitive comparison of the single 
DA maps is afforded by comparing their angular power spectra. In that case, it is easier to compare across frequencies 
because differences in beam response are readily accounted for by deconvolution. See Hinshaw et al. (2003a) for such 
a comparison. 

We generate three additional sets of difference maps using different combinations of the 4 channels within a DA. 
Specifically, we form the differences in the time-ordered data then generate maps as follows 

-(d!3+di4)--(d23+d24) A 12 , 

-(di3-di 4 )+-(d23“d24) “> A 34 , 

-(di3“di4)--(d23-d24) — ^ A 1234 , (27) 

where -4 indicates the map-making process. The A 12 maps are based on the polarization data, but processed as 
temperature maps, i.e., without attempting to demodulate the polarization signal. Since the two radiometers within 
a DA have completely independent detection chains, and since the polarization signal is weak to begin with (and is 
further suppressed by the lack of demodulation) the noise properties of the A 12 maps should be virtually identical to 
the nominal signal maps. The A 34 and A1234 maps are based on channel differences, (d /3 -d/ 4 ), and since the two 
channels within a radiometer have partially correlated noise, the noise properties of these latter maps will be different 
than the maps based on (d /3 + d/ 4 ). However, these maps do provide a check on the channel calibration, common-mode 
thermal effects and other potential artifacts. 

For each difference map the two-point correlation function and the angular power spectrum are calculated. The 
results are shown in Figures 10 and 1 1 and summarized in Table 6. Figure 10 shows the two-point function computed 
from the A J2 maps for Q2, V2, and W2. The most apparent feature in each of these functions is the slight bump at the 
beam separation angle of #beam ~ 141°, as expected (§2.1); the first data column of Table 6 gives C(0 beam )/C(O), for 
each DA. Note that, with the exception of K band, the ratio is typically 0.3%. The larger K band values arise because 
Galactic leakage in these difference maps is most severe in this band. This is also the source of the weak large-scale 
feature in the V2 two-point function in Figure 10. While this residual signal is small compared to the temperature 
signal, it is a systematic error that must be accounted for in the analysis of polarization data (Kogut et al. 2003). 

Figure 1 1 shows the angular power spectra of the difference maps for each of the 10 DAs, as well as for the final 
signal maps. Table 6 summarizes their statistics. Note that the white noise plateau in the 4 channel combinations per 
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DA divide into two families, as noted above, due to the correlations between channels 3 and 4. As a result, the null 
combinations, A34 and A 12 34, cannot be used to estimate the white noise parameter cr 0 for the signal maps. However, 
the polarization channel is seen to be in excellent agreement with the temperature channel in the white noise tail, thus, 
to the extent that real polarization signals, or other systematics, such as bandpass mismatch are not important, these 
maps should provide an excellent noise model for the temperature data. Table 10 summarizes the shape of the angular 
power spectrum at low l in the same way Table 5 did for the simulation. We find the spectra of these difference maps 
to be remarkably flat, with residual quadrupole moments of <4 /1K 2 for all bands except K (in which the difference 
is dominated by bandpass mismatch) and a single combination of W3. This value is much smaller than the small 
quadrupole measured in our sky (Bennett et al. 2003b). The power in the other multipole ranges is very close to 
the white noise floor, as seen in the final columns of Table 6. Since the residual signals seen in the flight difference 
maps are somewhat lower than those seen in the simulation (Table 5), we adopt the more conservative limits from the 
simulation as systematic error limits for structure at low /. This allows for the possibility than some of the error seen 
in the simulation comes from, e.g., common-mode calibration errors that cancel in the difference maps. 

The two-point correlation function of the A 12 maps demonstrates that the angle-averaged off-diagonal terms of 
the pixel-pixel covariance matrix are less than ~0.3%. However, the maps in Figure 6 indicate the potential for stripes 
along the scan directions for which the covariance can be locally higher than the angle- averaged value. In order 
to determine the magnitude of the covariance along the scan directions we perform the following computation. We 
form W band difference maps: Wi-Wi', where i - 1-4, and W i* is the average of the 3 other W band maps, e.g., 
Wl' = ^(W2+W3+W4). We then form time-ordered data from this map using the pointing appropriate to DA W i 
and compute the auto-correlation function, C(Af), from 30 days of data. This provides a measure of the pixel-pixel 
covariance along a stripe. The results for W3 and W4 are shown in Figure 12; the top panels show the covariance, 
normalized to one at lag zero, computed from the unfiltered maps, while the bottom panels show the results for the 
filtered maps. In order to convert the time axis to angular displacement along a scan, recall that the 2.°784 s _1 spin 
rate translates to a 2?6 s _1 rate for the beams in either focal plane (the second decimal place depends on position in 
the focal plane, and time in the precession cycle). The W3 result shows a hint of covariance (~0.1%) at lag 0.1 s, or 
0?26, prior to filtering, but none after (<0.05%). Prior to filtering, the W4 result shows clear covariance of up to 0.5% 
at small lag, decaying to <0.1% at lags of ~ 10 2 s, roughly one full spin. After filtering, the covariance is reduced 
by nearly a factor of two, but is still clearly detectable. This is understood to be residual covariance that survives the 
filtering process because of the fact that we subtract an estimated sky signal, based on the data, prior to filtering the 
noise, then add it back in to restore the signal. The algebra of this process is presented in Appendix B. We emphasize 
that W4 is the worst DA for 1 // stripes by at least a factor of 3 (Jarosik et al. 2003b), and we limit covariance along 
scan directions to be <0.1% for all other WMAP first-year sky maps. 


3.3. Pointing and Beam Determination 

3.3 A. Spacecraft Attitude Control and Determination 

The spacecraft attitude is determined from a combination of two autonomous star trackers (ASTs) with boresights 
perpendicular to the spin axis (along the spacecraft ±y axes), two rate gyroscopes, and two digital sun sensors. The 
sensor outputs are combined using a Kalman filter to determine the aspect solution. The sensor noise parameters and 
offsets were initially calibrated in flight during the in-orbit checkout (IOC) period in July 2001. By the end of IOC, 
the final tables were uploaded to the spacecraft. 

Spacecraft quaternions output by the Kalman filter provide the definitive transformation from the spacecraft 
reference frame to the J2000 geocentric inertial (GCI) system. Errors in the attitude solution are estimated using the 
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residuals of the individual sensor signals and propagated to the quaternions. After the final Kalman filter parameters 
were loaded, quaternion differences show a noise-like error with a 10" rras. In addition to the sensor noise, there 
is an apparent spin-synchronous error of ~ 10" that is believed to be due to propagation errors in the Kalman filter. 
As discussed below, this error is apparently seen in the instrument boresight determination using Jupiter observations. 
Since the pointing performance exceeds the requirement of 0!9 (root-sum-square for three axes), no correction of the 
spacecraft quaternions is attempted for the first-year processing. Sufficient information exists in the raw telemetry to 
attempt a correction in the future, if it is warranted. Note that random quaternion errors are automatically accounted 
for in the flight beam response maps generated from the Jupiter observations (Page et al. 2003a). 


33.2 . Instrument Boresight Determination 

As mentioned above, the spacecraft quaternions provide the definitive reference frame for the spacecraft. The 
instrument boresights, 10 each on the A and B sides, are determined from the Jupiter beam maps (Page et al. 2003a), 
which are generated with respect to the spacecraft frame provided by the quaternions. The boresight is defined as the 
location of the peak of a circular Gaussian fit to the main beam. The results of this fitting are given in Table 7 as 20 
unit vectors in spacecraft coordinates. These are the values used to determine instrument pointing in the first-year data 
processing. The uncertainty in the boresight position is ~ 2" per beam in both spacecraft azimuth and elevation. In 
addition to statistical uncertainty in the boresight fits, there are two other potential sources of error in the boresight 
determination: changes (drifts) with time, and errors in the relative time-tagging of quaternion data and science data. 

To test stability, we note that WMAP is in a position to see Jupiter twice per year for about 45 days each time. 
We refer to each ~45 day period as a Jupiter “season”. During the first season of each year, the boresights scan 
across Jupiter from roughly ecliptic north to ecliptic south, and vice-versa in the second season. As a test of boresight 
stability, we have generated beam maps from each of the first two seasons of data separately, and have fit boresight 
directions to each. We find the azimuth positions agree to better than 3" on both the A and B sides, but the elevation 
positions differ by ~ 10" on the A side, and a smaller amount on the B side. This difference between seasons is 
consistent with the ~ 10" spin-synchronous error in the spacecraft quaternions discussed above. We ignore this small 
effect in the first year processing, and subsume the small systematic error that results into our error estimate for the 
beam transfer functions, as discussed in Page et al. (2003a). 

The relative time-tag accuracy of telemetry packets was tested on the ground. A timing computer was set up to 
simultaneously receive test pulses from both the Attitude Control Electronics (ACE) box and the Digital Electronics 
Unit (DEU), the two computers that tag the attitude and science data packets, respectively. Each of these boxes in 
turn derives its time from the main “Mongoose” computer on WMAP (Bennett et al. 2003a). This test demonstrates 
a relative time-tag accuracy of 30 fi s between the quaternion packets and the science packets. In observing mode, the 
boresights sweep the sky at a rate of 2? 6 s" 1 , so a 30 fis time error produces a negligible pointing error of < 0.3" . 

Uncertainty in the spacecraft position is another potential source of boresight determination error. For the first 
year processing we use a predicted ephemeris that is uploaded to the spacecraft approximately weekly for on-board 
use by the Attitude Control System. This solution is returned in telemetry and is the basis for the ephemeris data 
supplied with the first-year release of time-ordered data. The uncertainty in these predictions are < 7 km in position 
and < 1 cm s~ ] in velocity (3 cr), relative to the Earth. An error of 7 km in WMAPs position would result in error of 
~2 mas in the apparent position of Jupiter as seen from WMAP and is thus completely negligible. 
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3.3.3. Beam and Window Function Determination 

Along with gain calibration and noise properties, knowledge of the beam shapes and window functions are among 
the most important aspects of the instrument to characterize for accurate measurements of the CMB. Page et al. (2003a) 
describe in detail the process by which beam maps are formed from in-flight observations of Jupiter, and how those 
maps are transformed to determine the beam window functions. The primary result they derive are a set of 10 beam 
transfer functions, bi , one per DA, based on azimuthally-averaged beam radial profiles. These transfer functions 
are included in the first-year data release (Limon et al. 2003). In addition, they derive a full covariance matrix for 
each transfer function, E^,, which characterizes the uncertainty in £>/. Typically, the uncertainty for a single DA is 
about 1-2%, with moderate covariance in /. See Figure 5 of Page et al. (2003a). These estimates already include a 
systematic error allowance to bound the small range of results obtained from different analysis methods. As described 
in Hinshaw et al. (2003a), the window function covariance matrices are propagated into the Fisher matrix (inverse 
covariance matrix) for the final combined angular power spectrum. Thus the final power spectrum, and the parameter 
fits based on it, already include statistical and systematic window function uncertainties (Spergel et al. 2003; Verde 
et al. 2003; Peiris et al. 2003). 


3.3.4. Effects from Elliptical Beams 

The WMAP beams are moderately elliptical, so the use of azimuthally-averaged radial profiles to describe the 
beam response is an approximation. This approximation is justified in Page et al. (2003a) by noting that the WMAP 
scan pattern produces excellent azimuthal averaging of the beam response in a large fraction of the sky. They have 
placed limits on the variation of the window function across the sky by comparing the effective window function in 
the ecliptic plane, based on a full two-dimensional transform of the beam response averaged over the flight range of 
scan angles, to the fully averaged transform, bi. For the three highest frequency cosmology bands, these variations 
range from 2-3% at Q band to ^1% in V and W bands. These variations are consistent with estimates of the angular 
power spectrum using data at high and low ecliptic latitudes separately (Hinshaw et al. 2003a). Since most of the 
statistical weight at high-Z resides in the V and W band data at high ecliptic latitudes, the use of fully averaged beam 
transforms is appropriate, and the systematic error estimate incorporated into should encompass any error in this 
approximation. 

Elliptical beams can also produce errors in the sky maps that are difficult to characterize in a simple way. We can 
define the sky map error due to non-circular beam response as 

Atasym = t 0 bs — tcirc, ( 28 ) 

where t 0 bs is the hypothetical noise-free sky map obtained with the actual experimental beam and scan pattern, and 
tcirc is the ideal sky map obtained by convolving the true sky with the averaged beam transform, bi . For a differential 
experiment like WMAP , there are two effects that contribute to At^m- The first, as noted, is incomplete azimuthal 
coverage in a given pixel, which gives rise to slightly elliptical peak structure at low ecliptic latitudes (see below), 
the second is a localized effect due to echoes from bright Galactic sources propagating to other pixels in the map. 
Specifically, as a bright source is observed in different orientations, the differential signal changes with orientation. 
Since the map-making algorithm must assign one average value to the pixel with the bright source, the ring of pair 
pixels at the beam separation will see an echo with a quadrupolar temperature distribution around the ring. We mitigate 
this effect by incorporating a bright source mask in the map-making algorithm, which is invoked as follows. If side A 
observes a pixel in the bright source mask, we only update the sky map accumulator for pixel A, but not for pixel B 


n<>bs * t n+ i (p A ) -> ^ • t„+i (p A ) + wi [d(t/) + 1 n (p B )] iW^) n^s (p A ) + w. 
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nobs • t n+ i(pB) A n ob ^t n+ i(p B )-w i [d(t i )-t n (p A )] n obs (p B ) A Robs (Pa)+vv t *, (29) 

where the terms are as defined in after equation (19). In this way we obtain an estimate of t (pa\ but we do not 
propagate bright echoes to the ring of neighbor pixels, of which ps is one. The mask we use for assigning this cut is 
the same Kp8 processing mask we used for the calibration fits (§2.2.1). 

We have generated a simulation to quantify the errors from both of these effects. Specifically, the simulation 
generates one year of noise-free differential sky signal which includes a model for the flight beam ellipticity. We run 
this data through the flight map-making pipeline to generate sky maps, t 0 b S . We also generate convolved maps t C ir C 
using the azimuthally-averaged beam transforms appropriate to the beam model. The residual map, At asym , for DA 
K1 is shown in Figure 13. The K band radiometers have the largest beam ellipticity of all the DAs, so this represents a 
worst case result. The general “mottling” near the ecliptic plane results from the relatively limited azimuthal coverage 
in this region producing elliptical peaks and anti-peaks which, in Atasym, are differenced with circular counterparts. 
This is especially noticeable near bright Galactic sources. The rms amplitude of these fluctuations in the Kp2 cut sky 
is 2 pK in K band, and at least a factor of 2 lower in Q-W bands. The effect of this structure in the power spectrum 
is primarily represented as a variation in the window function across the sky, as discussed above, and in Page et al. 
(2003a) and Hinshaw et al. (2003b). However, this structure also contributes to the 4-point function of the data in the 
sense that it couples power from different / ranges. This effect is potentially important for the interpretation of any 
gravitational lensing analysis of the WMAP data. 

The figure also exhibits faint echoes of the brightest sources that evade the map-making cut discussed above. We 
limit localized features in the Kp2 cut sky to less than 10 pK in K band and less than 2 pK in Q-W bands due to a 
combination of dimmer sources and more circular beams. We estimate that such features occupy <0.1% of the Kp2 
cut sky. 


3.3.5. Far Sidelobe Pick-up 

The WMAP optical system was designed to produce minimal pick-up from signals entering the far sidelobes. 
Barnes et al. (2003) present a complete determination of the WMAP sidelobe response by combining measurements 
from a variety of ground-based sources with in-flight measurements of the Moon. They produce response maps 
covering 47r sr that are included as part of the first-year data release. They then use these response maps, with the 
first-year sky maps, to estimate the systematic artifacts remaining in the first-year maps based on the well-justified 
assumption that sidelobe artifacts are small relative to the sky signal. The K band data have the largest sidelobe signal 
due both to the largest sidelobe spill, and to the brightest Galactic signal. The signal was deemed to be large enough, 
and well enough characterized, to warrant a small post-processing correction to the first-year K band map. Limits on 
remaining sidelobe induced artifacts in all the bands are presented in Table 1 of Barnes et al. (2003). 


3.4. Environmental Effects 

3.4. 1. Thermal Effects 

The radiometer gain and offset are dependent on temperature. There are several aspects of the WMAP design 
that are critical to mitigating this source of systematic error (Bennett et al. 2003a). The instrument is differential, so 
thermally induced gain changes act on a relatively small offset signal. The observatory environment was designed 
to be as stable as possible, consistent with other goals. For example, all nominal thermal control is passive to avoid 
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heaters cycling on and off. The observatory is placed at the second Earth-Sun Lagrange point far from the Earth, and 
the solar panels maintain a fixed 22?5 angle with respect to the Sun during normal observing mode. The instrument 
temperature is monitored with precision platinum resistance thermistors (PRTs) to verify the degree to which thermal 
stability is in fact achieved. 

Temperature variations at the spin period are the most critical since they can induce signals that couple relatively 
efficiently to the sky. However, owing to the relatively fast (129.3 s) spin period and the thermal mass of the instrument, 
any induced signals will have a very red spectrum and thus will couple only to the lowest few harmonic modes on 
the sky (/< 10). In the analysis below we use flight data to estimate the susceptibility of the gain and baseline to 
temperature variations of the instrument. In turn we use limits on the instrument’s physical temperature variation at 
the spin period from Jarosik et al. (2003b) to put limits on thermally induced artifacts in the time-ordered data, and 
hence the sky maps. 

Thermally induced signals can enter either through changes in the gain acting on the instrument offset or through 
changes in the offset itself. We show below that the latter are more significant for WMAP. 

The radiometer gain model presented by Jarosik et al. (2003b) describes the gain in terms of the RF bias (“total 
power”) of the detector, and the temperature of the FPA. This model tracks thermal variations in the gain on the time 
scale of the RF bias readout (23.04 s), and the map-making algorithm updates the gain on this time scale. However, 
since this sample rate is only a few times per spin, it is possible that a systematic temperature variation at the spin 
period could induce gain changes that are not well sampled by this model. As a separate check of gain-induced 
artifacts, we infer the temperature susceptibility of the gain from data taken over a long time period where gain 
changes are measurable. Results are given in Table 8. We combine these measurements with the limits on temperature 
modulation at the spin period derived by Jarosik et al. (2003b) to place limits on gain induced artifacts, as shown in 
Table 9. We conclude that gain-induced signals at the spin period, which might be poorly tracked by the gain model, 
are <20 nK. 

The instrument baseline is the product of the gain times the offset. As described in §2.2, we get an initial estimate 
of the baseline from the dipole calibration algorithm. This gives us an estimate of the instrument baseline once per 
hourly precession period. Sample hourly baselines for channels VI 13 and VI 14 are shown as a function of time over 
the first year of operation in Figure 14. Also shown is the temperature of the instrument FPA over the same time period; 
there is a clear temperature dependence in the baseline. We measure the baseline temperature susceptibility by fitting 
the hourly baseline estimates to a model of the form b(r) = co +c\* +C 2 ATfpa( 0 where the c ; are model coefficients and 
A7fpa = 7fpa - (7fpa) is the deviation of the FPA temperature from its mean. The most robust susceptibility results 
come from fitting a portion of the data near the time of a partial battery cell failure which occurred on day 2002:054 
(GMT) (Limon et al. 2003). In response to this event, the spacecraft bus voltage was autonomously commanded lower 
on day 2002:058 (GMT) causing the spacecraft to dissipate less power and thus cool slightly. The coefficients are given 
in Table 8. We have combined the results for the two channels in each radiometer because this is the combination that 
enters into the final sky maps. This has the effect of canceling some of the common-mode susceptibility measured in 
individual channels. As noted above, we combine these susceptibility measurements, taken over long time periods, 
with limits on the temperature variations at the spin period (Jarosik et al. 2003b) to place limits on induced signals at 
this time scale. The results are given in Table 9. We conclude that offset-induced signals at the spin period are <180 
nK. 

Slow drifts in the instrument temperature will be largely filtered out by the baseline pre-whitening discussed in 
§2.3.2. The steepest temperature gradient observed during the first year of observation occurred just after the above- 
mentioned battery cell failure. To assess the efficiency of the filtering process, we have analyzed the data during this 
period as follows. We applied the baseline pre- whitening filter to the temperature signal, 7 fpa( 0, to measure how much 
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the cooling gradient was suppressed by the filter. The input temperature gradient on day 2002:058 (GMT) was -1.7 
mK hr -1 . Applying the K1 1 filter to 7 fpa( 0 yielded an output gradient of -0.1 /x K hr -1 , while the W41 filter yielded 
an upper limit 1 0 times smaller. We conservatively estimate upper limits on residual drift in the filtered baseline of 
< 10 nK hr” 1 for the most susceptible channels. 


3.4.2. Electrical Effects 

A variable electrical signal on board the observatory could induce an apparent signal in the radiometers. Sources 
of variable electrical signals include the reaction wheels, transponder, bus voltage fluctuations, and RF noise coupled 
to the instrument through the power bus. During the final observatory thermal vacuum test, in which the observatory 
was operating at temperatures close to those achieved in flight, searches for such electrically induced radiometric 
artifacts were conducted (Jarosik et al. 2003a). Upper limits on radiometer bus voltage susceptibility, based on ground 
tests, are given in Table 8. We combine these results with an upper limit on bus voltage variations of 3.0 mV rms, 
measured on-orbit, to conclude that electrically-induced signals at the spin period are <40 nK. See Table 9. 


3.5. Miscellaneous Effects 

3.5.1. Radiometer Cross-talk 

A large signal in one radiometer could induce an erroneous signal at the output of another radiometer due to 
electrical cross talk. Such cross talk is not expected but could arise from, e.g., non-ideal amplifier behavior, or other 
parasitic effects, such as pickup in the wiring harnesses. 

A careful search was made for this effect during the instrument ground tests. Noise diodes were used to inject 
a large signal into one radiometer at a time while the input feeds of all other radiometers were covered by absorptive 
loads. The outputs of the 9 non-driven differencing assemblies were searched for any evidence of the injected signal. 
The tests were run with the amplifiers in the passive channels both on and off in order to distinguish pickup mecha- 
nisms. No pickup was found in any test. Table 10 gives 2 o upper limits on the pickup by any DA due to any of the 9 
other DAs. The column labeled Electrical gives the results obtained from the test with the amplifiers turned off, and 
the column labeled Radiometric gives weaker limits from the test with the HEMTs turned on. The latter limits are 
weaker because the output noise levels are higher. 

This level of pickup is quite small and could only be of potential concern when WMAP scans across Jupiter, the 
brightest source in the sky for WMAP at L2. The values in the Table 10 are thus referred to peak Jupiter signals in 
each band. For example, the first entry indicates that when Jupiter induces a signal of 185 mK in the W1 radiometer, 
the pickup in the K1 differencing assembly is <30 /xK (95% confidence) which is -26.8 dB below the peak Jupiter 
signal of 14 mK in K band. This signal occurs when the beam of the pickup channel is within a few degrees of Jupiter, 
depending on channel separation in the focal plane, and in every case is less than direct radiometric detection of Jupiter 
in the near side lobes. 

This limit on cross talk implies pickup is completely inconsequential in normal observing mode. Using the same 
Table example, a 100 /xK signal in W1 could cause at most a 200 nK signal in the most susceptible the four K band 
differencing assemblies. 
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3.5.2. Source Variability 

Time variable objects are a potential source of contamination for observations of the CMB; see for example 
Sokasian et al. (2001) and references therein. One concern is that an object may grow in brightness over the course of 
WMAP's observations, avoid detection during an initial source survey, and remain unmasked during subsequent data 
analysis. For example, blazars produce relatively rapid and large amplitude variability in all wavebands. Long term 
observations of such objects show that increases in flux by a factor of up to a few over a time scale of years can be 
anticipated (Flett & Henderson 1983; Ennis et al. 1982; Stevens et al. 1994; Bower et al. 1997). Observations of Zw 
2 by Falcke et al. (1999), provide an extreme example: a greater than 20-fold increase in brightness, from ~0.1 Jy to 
~2 Jy, over a period of less than two years. While this object could produce a temperature response of a few hundred 
/xK in the WMAP data, such objects are rare and, if left undetected, would have a minimal effect on cosmological 
inferences. Tests for point source contamination in the WMAP data are given by Bennett et al. (2003c) and Hinshaw 
et al. (2003b). These tests will need to be revisited on an annual basis. 

Another source of concern is that a time variable source in the WMAP data has the effect of broadcasting noise 
to the ring of ~1000 pixels which are separated from the variable source by the beam separation angle (0beam ~ 141°). 
The point source list derived from the WMAP first-year data is 98% reliable with ~5 spurious detections at the ~0.5 
Jy flux limit of the survey (Bennett et al. 2003c). The nominal point source sensitivity of the WMAP telescope is 
T ~ 200 /xK/Jy, thus a noise level of ~0.1 jjK is expected from variable sources that evade detection. 


4. CONCLUSIONS 

The processing steps used to produce the first-year WMAP sky maps include an initial simultaneous estimate of 
the sky map and the instrument calibration. The instrument gain is then refined using a model based on engineering 
telemetry, and the instrument baseline is refined by the application of a pre- whitening filter. A final archive of calibrated 
data is produced and used to generate final sky maps using a slightly refined iterative algorithm. 

We limit systematic artifacts due to calibration, map-making and environmental disturbances to less than ~15 
fjK 2 in the quadrupole C 2 , with tighter limits at higher multipole moments (Table 5). Beam transfer functions are 
measured for each beam with 1-3% over the entire range of multipole moments that WMAP is sensitive to (Page et al. 
2003a). The covariance matrix of the beam transfer function is propagated through to the final power spectrum error 
matrix. We characterize pixel-pixel covariance matrix and place limits on residual stripes in the final maps. 

All major data products from the first year of WMAP observations are being released through NASA’s new 
Legacy Archive for Microwave Background Data Analysis (LAMBDA) at http : / / lambda . gs f c . nasa . gov/ . 


The WMAP mission is made possible by the support of the Office of Space Sciences at NASA Headquarters and 
by the hard and capable work of scores of scientists, engineers, technicians, machinists, data analysts, budget analysts, 
managers, administrative staff, and reviewers. We acknowledge use of the HEALPix package. 
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A. THE MAPPING FUNCTION 


Equation (4) defines the continuous form of the mapping function, which encodes both the scan strategy of an 
experiment, and convolution due to the beam response. We can relate this to the matrix form, in equation (7) as 
follows. The mapping function evaluated at time for a finite integration time r may be written in terms of the beam 
response function as 

M(n,f,)= - [‘ T dt [a (R(r ) ■ n) - /? £#(R(r) • n)] , (Al) 

T J<> 

where B^(n) is the beam response of the A-side beam, in spacecraft coordinates, normalized to unit integral 


Jdn„B A (n)=l, (A2) 

similarly for the B side, and R(r) is the rotation matrix from sky-fixed (Galactic) coordinates to spacecraft coordinates 
at time t . The terms a and p in equation (Al) account for possible ohmic losses in the A and B-side optics that are 
not necessarily equal (Jarosik et al. 2003b). Since the data are calibrated using the modulation of the CMB dipole, we 
only need to parameterize the loss imbalance which, following Jarosik et al. (2003b), we parameterize as 


a = 1 +* im , 

P = l-x im . (A3) 

Note that loss imbalance is separate from lossless differences in the beam response function, e.g. differences in the 
solid angle of the A and B-side beams. Once the calibration is applied, the differential sky signal is a measurement of 
the form 

At(t) = Jd^ n t(n)[(l+x im )B A (R(t)-n)-(l-x im )B B (R(t)-n)), (A4) 

which still includes the effects of any loss imbalance. We now separately consider how this calibrated differential data 
propagates into the sky maps and the Jupiter beam maps. 

When making sky maps from the calibrated data, each datum is modeled simply as 

At(f,) = (1 +*im)t(/? / i)-(l -X im )t(p B ), (A5) 

where is the pixel observed by the A-side beam at time t if and similarly for ps. That is, each row of the mapping 

matrix in equation (7) has the form 


M(/?,f|) — (1 +Xi m ), ... , (1 Xi m ),...], 


(A6) 


This preprint was prepared with the AAS macros v5.0. 
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with non-zero entries in pixel columns p A and p B only. Upon solving for the sky map, this ideally leads to an effective 
beam response of the form 



_fi£V) 1 flgV) ! 

(l+*im) (1-^im) 


(A7) 


where is the symmetrized beam response for the A-side beam and similarly for the B-side. For this ideal case 
to obtain, the following must hold: 1) each pixel must be observed equally by the A and B-side beams, which is 
true to a very good approximation for WMAP , and 2) each pixel must be observed with uniform azimuthal coverage. 
Deviations from these assumptions are considered in the text. 


The beam mapping data is compiled from calibrated observations of the bright source Jupiter. The calibrated data 
taken when side A is observing Jupiter has the form 


At — J <m n tjr(n) ( 1 +x im ) B A (R*n)+ Ats ky , (A8) 

where tj(n) is the brightness temperature of Jupiter in the direction n and At sky is the background sky temperature 
difference, which is subtracted during processing. An analogous equation holds for the B-side data. Assuming the 
beam response is constant over the extent of Jupiter, the integral reduces to 


At — Tj Qij ( 1 + Xim) B a (R • n/) 5 + At sky , 


(A9) 


where 7} is the disk brightness temperature of Jupiter, and flj is its solid angle. Beam maps are compiled by binning 
the corrected data At- At sky as a function of n, in spacecraft coordinates. This produces maps proportional to the 
beam response 

T A (n) = TjnXl+Xim)B A (n), (A 10) 

and similarly for the B-side. 

Ultimately, we wish to compute the transfer function of the symmetrized beam response. This may be obtained 
from the symmetrized beam maps as 


B(n) = 


1 

2 Tjtlj 


7f(n) | 7f(n) 

( 1 + Xj m ) ( 1 — *im) 


(All) 


We don’t know the brightness temperature of Jupiter a priori , but since this is an overall normalization factor, we are 
free to normalize the final transfer function to 1 at / = 0. 


B. MAP-MAKING WITH FILTERED DATA 

§2.3.2 presents the filtering algorithm used to determine the final instrument baseline. This process included an 
estimated sky signal subtraction based on the initial sky maps produced with the hourly calibration. In this Appendix 
we derive the noise properties of sky maps produced with this filtered data. In the following, we assume the time- 
ordered data has a noise covariance N = (nn r ) that includes a 1 // component, and that we have a pre- whitening filter 
F such that 


<(Fn)(Fn) r )ocI. 


(Bl) 
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B.l. Map-making with filtered signal + noise 

We could filter the full data prior to any calibration or sky map estimation, then deconvolve the effects of the filter 
in the subsequent data processing. The input data would have the form 

d' = Fd = FMt+Fn. (B2) 

Then, in order to obtain an unbiased sky map estimate, we would need to evaluate the sky map estimator 

t' = (M 7 'F r FM)" 1 • (M r F r d'), (B3) 

which deconvolves the action of the filter on the sky signal. Since the 1 // noise in the WMAPdata is relatively small, 
implementing this estimator for the first-year sky maps was deemed unnecessary and would likely have delayed the 
release of the maps. The alternative is to filter only the noise by subtracting an estimate of the sky signal prior to 
filtering, then adding it back in to the time-ordered data prior to making new maps. 


B.2. Map-making with filtered noise 

Let to be the sky map estimated from unfiltered data, using the hourly calibration. This is related to the true sky 
signal by 


Wd 

04) 

W-(Mt+n) 

(B5) 

t+Wn 

06) 


where W = (M^M)" 1 * M T is the map-making operator defined in §2. 1 , and we have used the fact that W ■ M = 1 . We 
use this sky map to subtract a sky signal from the time-ordered data prior to filtering, then we add it back in after 
filtering. This produces a filtered data set 


F • (d- Mto)+Mto 

(B7) 

F • (Mt + n - Mt- MWn) + Mt+MWn 

08) 

Mt+Fn+(1-F) - MWn. 

(B9) 


This time series data set consists of an unbiased sky signal Mt, a white noise term Fn, and a residual noise term 
(1-F) ■ MWn which is due to the off-diagonal “wings” of the filter (1-F) acting on the noise from the initial sky map 
estimate, MWn. 

We can make a map ti from the data di using the algorithm of §2.1 


Wd, 

(BIO) 

t+WFn+W (1-F) MWn 

(Bll) 

t+WFn+(l-WFM)Wn 

012) 

t+WFn+R Wn, 

(B13) 


where we have again used W • M = 1 and we have defined a “residual” operator R = (1 - WFM) which is small in the 
sense that only off-diagonal terms in F contribute to it. This is most easily seen if we recall that the filter operator is 
1 on the diagonal and has small off-diagonal terms. We can then write F = 1-E from which it follows R = WEM. 


I 
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It follows that ti is an unbiased estimate of t that includes a white noise term WFn (this noise still contains the small 
beam separation covariance), and a residual noise term, R * Wn, due to the noisy sky signal estimator used in the 
filtering process. This latter term is the “excess” noise seen in the W band single DA maps after filtering (Hinshaw 
et al. 2003a). 

The residual noise term can be reduced somewhat by iterating the filter algorithm a second time, using 1 1 as a sky 
signal estimator. It is straightforward to show that, after some algebra, the resulting map is 

t 2 = t+WFn+R-WFn+0(R 2 ), (B14) 

where we have neglected a term of order R 2 • Wn. The residual noise is reduced slightly since R is now acting 
on the white noise WFn instead of the full noise Wn. Rut since the white noise dominates, this is a relatively 
insignificant improvement. It is also clear that subsequent iterations of the filtering only contribute higher order 
corrections which are negligible. This convergence has been verified with the flight data. Another feature we have 
verified with the flight data is the fact that the excess noise term decreases with time. The reason for this is simply 
that R = (1-WFM) = WEM gets smaller with additional data because the map-making operator W gets smaller as 
more observations accumulate in its “denominator”, (M 7 M) _1 . The first-year maps were only filtered once, as per 
equation (B13), because the improvement in noise properties was not deemed sufficient to warrant the additional 
processing time. 


C. MAP-MAKING WITH POLARIZATION 

WMAP observes the sky with two orthogonal linear polarization modes per feed, thus it is sensitive to the 3 
Stokes parameters 7, Q , and U. This Appendix outlines the algorithm with which these parameters can be mapped 
using the differential data from WMAP : The approach is an extension of the iterative method in §2.1.1 introduced by 
Wright etal. (1996). 


C.l. Polarization Mapping with Total Power Data 

Suppose we observed the sky with a single beam, total power radiometer that is sensitive to a single linear 
polarization, denoted mode #1. In terms of the Stokes parameters, the temperature observed by the instrument at time 
t in pixel p would be 

di(t) = i(p)+q(p)cos2'y+u(p)sin2y, (Cl) 

where i, q, and u are Stokes parameter maps in units of temperature, and 7 is the angle between the polarization axis 
of the beam and the chosen reference direction for pixel p (the choice of reference direction is discussed below). Note 
that we adopt the common convention that the instrument response reduces to the total temperature in the limit of an 
unpolarized source. 

The signal in the orthogonal polarization channel, which is fed by the other port of the ortho-mode transducer 
(OMT) and denoted mode #2, is given by 


d 2 (?) = i(/?) + q(p) cos 2(7 + ^ ) + n(p) sin 2(7 + ^ ) 
= i(p) - q(p) coslj -u(p) sin 27 . 


(C2) 

(C3) 
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We can isolate the intensity and polarization signals by taking sums and differences 

d(0 s ^(d,+d 2 ) = i(p) 

P(0 = | (d i - d 2 ) = q(p) cos 2 7 + u (p) sin 


(C4) 

(C5) 


Given noisy data, we can estimate the intensity i(p) by averaging all the data dfo). For the polarization we can only 
estimate a linear combination of q and u from a single observation. However, if we have several observations of pixel 
p with a variety of polarization angles 7 , we can estimate q and u in a given pixel by minimizing x 2 > defined as 


2 — 


=£ 


[pfe) ~ q(p) cos 27/ - u (p) sin 27/f 


(C 6 ) 


i€p 


where i is a time-ordered data index and the sum is over observations within pixel p, and 7 , is the polarization angle 
for the I th observation. The best-fit values for q and u are given by 


( q(/>) \ = ]_f -E,<W CT < 2 

V u (p) ) A V ~E t CiSi/af Ei C J l a l 


\ ( \ 

) V sE^/Pfc )/^ 2 ) ' 


(C7) 


where a = cos 27 /, $/ = sin 27 /, and A = c? /of Yli ~ C /of ) 2 is the determinant of the normal equations 
matrix. The standard errors for q and u are given by the inverse of the normal equations matrix 



/ 


(C 8 ) 

(C9) 

(CIO) 


In the limit of uniform azimuthal coverage and constant noise per observation (a; = ao), the matrix elements in the 
linear system reduce to 


( Ei c ? E/*m \,{ N / 2 0 \ 

V E/^, E ,*} ) V 0 N/2 ) ’ 


(Cll) 


where TV is the number of observations of pixel p. In this limit, the noise in q and u is equal and uncorrelated and 
reduces to 



(Cl 2) 


Thus the noise in each polarization component is \fl times noisier than in the intensity i. 


C.2. Polarization Mapping with Differential Data 

We now generalize to the case of polarization mapping with differential input data. For clarity, we first consider 
the case where the loss in the two sides A and B are equal. We generalize to the case with unbalanced loss in the next 
subsection. In the case of WMAP, one radiometer in a differencing assembly (radiometer #1) is fed from the axial 
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port of the OMT and the other (#2) is fed by the lateral port (Jarosik et al. 2003a). Following 
differential signals from radiometers 1 and 2 are 

equations (C1,C3) the 

di(r)= i(di3 + d] 4 ) = i(p A )+q(p A )cos2y A + u(p A )sm2'y A 

(Cl 3) 

- KPb) ~ <1(Pb) cos 275 - u(/?b) sin 2jb , 

(Cl 4) 

and 


d 2 (t)= i(d 2 3 + d 2 4) = 'KPa)-<\(Pa) cos 27 a - u(p A )sin 2^ A 

(Cl 5) 

- i(p B )+q(/7 B )cos27B + u(pB)sin27B, 

(Cl 6) 


where 7 a is the angle between the axial polarization plane and the reference direction in the pixel seen by the A beam, 
and similarly for 75. We take sums and differences of the two signals to isolate the unpolarized and polarized portions 
of the signal 

d(f) = i(di+d2) = ifo)-i(p») (C17) 

p(0 = ~(di-d 2 ) = q(/?A)cos27A + u(/?A)sin27A-q(p£>cos27 B -u(pB)sin27 B . . (C18) 


An iterative scheme for making maps of q and u follows the form used for intensity maps (§2.1.1). Let q„ and 
u n be the estimates of q and u respectively. Estimates for q n+ \ and u n+ i are obtained by combining the per-pixel 
fitting algorithm in equation (C6) with the iterative map-making algorithm, as follows 


* 2 =£ 

i€p 


[p'(fj) - q«+i (p) cos 27,- - u n+1 (p) sin 2jl \ 2 


(Cl 9) 


where the sum is over all observations of pixel p by either the A- or B-side beam, and p'(r) is the polarization data 
corrected with an estimate of the signal in the opposite beam 


f p(h) + QniPB ) cos 2 tb + u„ (pb) sin 27# beam A € p 

1 -p(t/)+qn(^A)cos27A +u 71 (/? A )sin27A beam B G p. 


The best-fit solution for q„ + i and u n+l is then 

fqn+i(p))_±f -H.c.Si/aj 

V u„ + ,(/>) ) A V -E iCiSi/cj 


Ei^p'fe)/®? ) 
E^/P'fc)/®? ) ’ 


(C20) 


(C21) 


where the sum on i is as defined above. The uncertainties are as given for the total power case, where <7/ is now the 
uncertainty per differential observation, p (*/). 


C.3. Map Making with Unbalanced Differential Data 

We now generalize to the case of map-making with unbalanced differential input data. In this case, losses in the 
A and B-side front ends are different and the differential signal is of the form 

(C22) 
(C23) 


d](0 = (l+^i m ,l)[i(/?A) + S(/?A ) 7A)]~(l-^im,l)[i(/?B) + S(/7 B? 7B)] 

d 2 (r) = (1 +Xim,2)[i(PA)“S(pA,7A)]~(l , 
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where *i m ,i,Xi mj 2 are the loss imbalance parameters in radiometers 1 and 2, respectively, as defined in Jarosik et al. 
(2003b), and s(p,7) is short-hand for the linear combination of Stokes parameters 

s(p,7) = q(/?)cos27 + u(p)sin27. (C24) 

As before, we take sums and differences of the two signals to isolate the unpolarized and polarized portions of the 
signal. First, define the mean imbalance and the “imbalance in the imbalance” as 

2 C^im, 1 + -*1111,2) (C25) 

2(*im,l (C26) 

then 


*im = 

Sx im = 


d(/)=l(d,+d 2 ) 

P(f)=^(d,-d 2 ) 


KPa) - i (pb) +*i in [i(/?a) + i (pb)] + Sxim [s (pa ,7 a)+ S (pb, 7b)] 

(1 +-*im)i(PA) — 0 ”^im)KpB) + CK^*im) 

S (pA ,7a)“ S (PB , 7 b) + *im [S (pa ,7b) + S (pb, 7a)] + Sx im [i (pa) + i(pa)l 
(1 +^m)s(p^,7 i 4)-(l -^im)s(pa,7B) + 0(Sx im ). 


(C27) 

(C28) 


Note that the term of D(Sx im ) in equation (C27) is negligible, because i » s, but the term in equation (C28) must be 
considered more carefully. First note that &ti m is small - from Table 3 in Jarosik et al. (2003b), the largest value is 
0.35% in W2, with most values being ~0.1%. Second, the multiplier, (i(/?A) + i(pa)) does not modulate with polariza- 
tion angle, 7, thus it is effectively an offset term that is further suppressed by the map-making demodulation. Finally, 
the term is out of phase with the differential signal (Kpa)~KPb)) so it does not effectively couple to the sky. The effect 
of this term in the first-year data is further considered by Kogut et al. (2003). 

We generalize the differential map-making algorithm to account for loss imbalance as follows. For intensity, let 
i n be the n th estimate of i, then 


f [+d(^) + (l-^im)in(PB)]/(l+^im) beam A £ p 
1 f-dft) + ( 1 + -Xim) inipA )]/(!“ Xim) beam B £ p. 


The updated intensity map is then estimated by binning the corrected data 


Wp) = 


£fW,-d # fo) 

E-'W,- ’ 


where w, is the normalized statistical weight of each observation 


f (l+*im) Vo/of beam A £ p 
1 (l~*im ) 2 OqI<t} beamB £ p. 


(C29) 


(C30) 


(C31) 


For polarization, let q„ and u„ be the rc th estimates of q and u respectively, and define s n (p,7) = q„(/?)cos27 + 
u„(p)sin27. Estimates for q„ +] and u„+i are obtained by combining the per-pixel q and u demodulation with the 
iterative map-making algorithm. We define \ 2 as follows 

[p'fo) - <Wi (p) cos 27/ - u„+i (p) sin 27/] 2 



(C32) 
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where the sum is over all observations of pixel p by either the A- or B-side beam, and 


/f , = f [+P(fi)+(l-^im)Sn(PB,7s)]/(l+^m) beam A £ p 
~\ [-p(r,-) + ( 1 + ^™)s n (pA , 7 a )] /( 1 -*im) beam B € p. 


The best-fit solution for q„ +l and u n+i is then 


( q n+ i (p) \ = W 

' I >*? 


( E /W/C, 

ip'fo) A 

V u n+ i(p) J A \ 

„ -E/WiCtfi 

E M ) 

V E ,-WiSi 

p'(fi ) ) 


where the weight w/ is the same as equation (C31). 


(C33) 


(C34) 


C.4. The Choice of Reference Direction 


We adopt the sign convention of Zaldarriaga & Seljak ( 1 997) in which the polarization components are defined in 
a right-handed coordinate system with the z-axis pointed outward towards the sky. The Stokes parameters are defined 
with respect to a fiducial direction in each pixel on the sky. We adopt the convention in which the reference direction 
is aligned with the local Galactic meridian, i.e., the great circle connecting a given point to the North Galactic Pole. 
The unit vector, n, tangent to this great circle, pointing in the direction of the North Pole, is given by 


n=l xe=lx 


(z x i) 
sin 0 


where 1 is the unit vector along the line of sight of the current observation, e is a unit vector pointing east from /, z is 
the unit vector to the North Galactic Pole, and 0 is the polar angle (co-latitude) of L 

For reference, we give formulae for computing the factors cos 27 and sin 27 here. Let / be the unit vector along 
the line of sight, w be the unit vector pointing west from /, h be the unit vector pointing north from l (the polarization 
reference direction) and p be the unit vector along the polarization plane defined by the axial port of the OMT. Then, 
for both the A and B sides, we have 


C0S7 
sin 7 
cos 27 
sin 27 


p'W 

2cos 2 7-l 

2sin7COS7. 


Note that this defines a right-handed coordinate system with (x,y,z) axes (n, vv, /) whose z axis is oriented outward 
following the sign conventions of Zaldarriaga & Seljak (1997). 
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Fig. 1 . — A residual sky map, t out - t in , from an “ideal” one-year simulation of Q2 data, designed to test the iterative 
map-making algorithm presented in §2.1. The input sky map included realistic CMB signal with a peak-to-peak 
amplitude of ~ ±420 /z K, and a Galactic signal with a peak brightness of ~50 mK. The rms structure in this map is 
<0.2 /zK, after accounting for the 0.15 /zK noise that was introduced to the simulation to dither the digitized signal. 
The map is projected in ecliptic coordinates and shows the anisotropy mode that is least well measured by WMAP, 
due to a combination of the scan pattern and the beam separation angle. This residual level is the result of 50 iterations 
of the algorithm - more iterations would reduce it even further. 
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Fig. 2. — Convergence of the dipole-based gain solution for a selected WMAP radiometer channel (K1 13) based on a 
one-year simulation. This simulation was generated with an input gain of 300.0 du mK“ J , and minimal noise. The first 
iteration, which assumes the sky model has only a dipole component, leaves residual gain errors of up to 7%, due to 
the projection of the relatively bright Galactic emission onto the dipole model. After 30 iterations of the simultaneous 
fit described in §2.2.1, the residual errors in the gain solution are less than 0.1%. 
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Fig. 3. — The hourly gain and baseline fit described in §2.2.1 from the flight data for channels K1 13 and VI 13. The 
top two panels show the gain solution, the bottom two the baseline. Note that the gain is stable to ~5% over the first 
year (see also Table 3). The variable noise is due to the changing projection of the scan pattern on the CMB dipole 
over the course of a year. The instrument baselines have a typical drift of 5-10 mK over the year. The channel VI 13 
exhibits one of the clearest thermal susceptibilities of the WMAP radiometers, though we show in §3.4.1 that the 
induced systematic signal is negligible. See also Figure 14. 
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Fig. 4. — The left-hand panels show the measured auto-correlation function, C(Ar)/C(0), for selected radiometers, 
of WMAP time-ordered data, after subtracting a model sky signal based on the initial sky maps. The model fits are 
indicated by an x at a lag of 1 observation and by straight lines for A t>r (see §2.3.2). All of the WMAPDAs except 
W4 have a covariance of < 1% at non-zero lag (Table 4). The anti-correlation at lag ~2000 s is due to the subtraction 
of the hourly baseline as a pre-filter. The middle panels show the pre-whitening filter, in the frequency domain, that 
is applied to the time-ordered data after a model sky signal has been subtracted. The vertical dashed line indicates the 
spin frequency, 7.7 mHz, and the number indicates the fraction of power transmitted by the filter at the spin frequency. 
The right-hand panels show the measured auto-correlation function for selected channels of WMAP time-ordered data, 
after pre-whitening, on the same scale as the left panels. The apparent change in noise level at different lags in C( Ar) 
is due to a step-wise change of bin size in At. 
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Fig. 5. — The time-ordered data for channel W424 before (black) and after (red) applying the pre- whitening filter. 
The data are boxcar averaged over a 46.08 s window to show the low frequency noise in the unfiltered data. Without 
averaging, the data before and after filtering are virtually indistinguishable. Note that baseline variations in this channel 
are of order 2 mK on a one-hour time scale, as expected given the measured 1 // knee frequency of this radiometer 
(Jarosik et al. 2003b). W4 is the worst differencing assembly from the standpoint of 1 // noise. 
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Fig. 6. — Simulated and flight difference maps showing the structure that is removed from the maps by the pre- 
whitening filter. All four maps are differences between sky maps generated before and after baseline filtering. The 
maps are projected in Galactic coordinates and the temperature scale on each is ±50 /zK. The “blobs” of white noise 
along the ecliptic plane can be ignored. They arise from differences in the handling of planet flags in the two forms of 
the map-making code. The top two panels show W3 and W4 data from a one-year simulation that includes flight-like 
1// noise in the time-ordered data. The bottom two panels are the same for the flight W3 and W4 maps. Note the very 
different structure between W3 and W4, due to different 1 // knee frequencies (Jarosik et al. 2003b). Note also that 
the simulation captures the basic structure of the flight data very well. 
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Fig. 7. — The recovered gain solutions for channel V 1 1 3 in a flight-like simulation (top) and in the flight data (bottom). 
The “noisy” black traces show the hourly baseline binned in 24-hr samples (to reduce noise) and the green traces are 
the best fit gain model (§2.3.1). For the simulation, the input gain used to generate the simulated data is shown in grey. 
In the simulation, the absolute gain is recovered to better than 0.1% in all 40 channels, and the binned hourly gain is 
everywhere within ~0.2% of the gain model, and the input gain. Gain changes are well tracked by the pipeline. 



-40- 



Fig. 8. — Angular power spectra of the 10 resdidual maps t out -t in generated from the flight-like one-year simulation. 
In each case the spectra were evaluated in the Kp2 cut sky (Bennett et al. 2003c). Table 5 quantifies structure in these 
maps beyond flat white noise. Note that most features are restricted to / * 25 but with an amplitude that is still much 
less than the sky signal in this range. The residual effects of 1 // noise are seen in the gradual rise of the noise spectrum 
at low / in W4. See Hinshaw et al. (2003b) for further discussion of this. 



-41- 



Fig. 9. — Sum and difference maps generated from the flight Q, V, and W band data, as indicated. To reduce the noise, 
all maps have been binned in larger pixels (HEALPix N S ide = 64) and displayed with a temperature scale of ±100 /zK. 
As discussed in the §3.2, the only apparent structure in the difference maps is due to residual Galactic contamination 
owing to the fact that the effective frequencies of the DAs are slightly different. This does not affect signals with a 
CMB spectrum, because the calibration source (the CMB dipole) has the same spectrum. See Bennett et al. (2003b) 
for higher resolution images of the signal maps. 
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Fig. 1 1. — Angular power spectra of signal and noise maps for each DA. In each panel, the upper red and green traces 
are the spectra of the null maps, A 34 and A 1234, respectively. The lower green trace is the A 12 map, and the black 
trace is the signal map. The blue curve is our best estimate of the underlying CMB signal from Hinshaw et al. (2003b). 
The pairing of white noise levels is discussed in §3.2, Table 6 presents a measure of structure in the difference spectra, 
which are remarkably fiat. 













Fig. 13. — Residual map from a K band elliptical beam simulation. The output map was generated from a one- 
year simulation of data with an elliptical beam response. The residual map shown was generated by subtracting the 
underlying sky signal convolved with the nearest effective circular beam response. This remaining structure contributes 
to the four-point fluctuation spectrum. The scale of the color range is ± 10 jj,K. The rms structure in the Kp2 cut sky is 
2 /iK. See §3.3.4. 
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Fig. 14. — An illustration of the mild thermal susceptibility of the instrument baseline. The top panel shows the tem- 
perature of the instrument Focal Plane Assembly (FPA) over the course of the first year. The second and third panels 
show the hourly baseline solution for channels VI 1 3 and V 1 1 4, which are among the most thermally susceptible. Note 
that the thermal baseline response is mostly common-mode. The channel combination that contains sky signal is the 
difference between channels 3 and 4, thus most of this response cancels. On day 2002:054 (GMT) a partial battery 
cell failure led to a commanded decrease in spacecraft bus voltage with a corresponding decrease in overall power 
dissipation and spacecraft temperature. This event provides a clean measurement of the instrument baseline thermal 
susceptibility - the bottom panel shows a close-up of the VI 14 baseline near this event. The dashed line is a fit to 
a model including a term proportional to db/dTppA . The best-fit susceptibility values for all channels are given in 
Table 8. See Limon et al. (2003) for a complete discussion of WMAPs first-year thermal profile. 
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Table 1 . WMAP data processing notation 


Symbol(s) 


Description 


Uh 

r 

N P 

N, 

P 

Pa,Pb 

$beam 

7 

7a. 7 b 

C|> Si, CiSi 

t (P) 

Kp) 

q (p), u (p) 

t ip), tn(p) 

nobs ip) 
tc (p) 
tg(p) 
*s(P,0 

At(r) 

Atd(r) 

Atv(r) 

Ata(r) 

c(0 

dy(0 

d(f) 

p(f) 

g(0 

n(0 

b(0 
g*. h 
o< 


(TO 

•^im» x imJ 

C(At) 

Cx,A,B 

w(f) 


Time, time of ith observation, in s 
Integration time per observation, in s 
Number of pixels in a map, 0 - (N p - 1 ) 

Number of time-ordered data points 
HEALPix pixel number 
A & B-side pixels at time t 
Separation of A & B-side beams, in degrees 
Polarization angle 

A & B-side polarization angles, at time t 
cos 27, sin 27, cos 27 • sin 27, at time u 
Sky map, in mK 

Sky map, Stokes parameter /, in mK 
Sky map, Stokes parameters Q , U , in mK 
Estimated sky maps, in mK 
Number of observation of pixel p 
CMB anisotropy map, in mK 
Galactic foreground map, in mK 
Time-dependent source map (Sun, Earth, Moon), in mK 
Time-ordered differential sky signal, in mK 
Time-ordered COBE dipole signal, in mK 
Time-ordered local velocity dipole signal, in mK 
Time-ordered anisotropy signal, At- Atd, in mK 
Time-ordered raw data, single channel, in du 
Time-ordered data, radiometer i, channel 7, in mK 
Time-ordered intensity data, co-added channels, in mK 
Time-ordered polarization data, co-added channels, in mK 
Rad. responsivity (oc gain), single channel, in du mK" 1 
Rad. noise, single channel in calibration, 

4 co-added channels in map-making, in mK 

Rad. baseline, single channel, in du 

Hourly gain, baseline, single channel, k^precession 

rms noise, /^observation, single channel in calibration, 

4 co-added channels in map-making, in mK 
Mean rms noise, single or co-added channels, in mK 
Loss imbalance parameter (radiometer j) 

Auto-correlation function of noise, in mK 2 
Auto-correlation function model parameters 
Pre-whitening filter, Fourier space 



Table 1 — Continued 


Symbol(s) Description 

F Pre-whitening filter, time domain, N t x N { matrix 

N Time-ordered noise covariance, N t x N t matrix, in mK 2 

M Mapping function, N t x N p matrix 

W Map-making operator, (M r M)~ l • M r ,N p xN t matrix 

S Pixel-pixel noise covariance, N p x N p matrix, in mK 2 

D Reduced inverse noise, (M T M) = Oq S _1 ,N p xN p matrix 

A ij Difference map from channel combination i j 
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Table 2. WMAP Attitude Control System Requirements 


Parameter 

Requirement 

Performance 

Precession rate (< d<f>/dt ) 

-0? Is -1 ±6.3% 

-0? Is' 1 ±3.6% 

Spin rate (dip/dt) 

2?784s-'±5% 

2P784s _1 ±0.13% 

Sun-spin angle (0) 

22?5±0?25 

22?5 ±0?023 


Table 3. WMAP Dipole-Based Gain Summary 


Channel 3 

Radiometer (g k ) A g* a 

(dumK" 1 ) (%) 

(%) 

Channel 4 
(gk) Ag/ 

(dumK’ 1 ) (%) 

(%) 

Kll 

-0.974 

7.4 

0.66 

+0.997 

6.8 

0.66 

K12 

+1.177 

6.2 

0.75 

-1.122 

6.4 

0.75 

Kail 

+0.849 

4.7 

0.75 

-0.858 

5.0 

0.75 

Kal2 

-1.071 

5.1 

0.75 

+0.985 

5.3 

0.75 

Qll 

+1.015 

4.5 

0.94 

-0.948 

4.4 

0.94 

Q12 

+0.475 

5.3 

1.03 

-0.518 

5.5 

1.03 

Q21 

-0.958 

5.8 

0.94 

+0.986 

6.0 

0.94 

Q22 

-0.783 

2.9 

1.22 

+0.760 

2.8 

1.22 

Vll 

+0.449 

4.5 

1.50 

-0.494 

4.5 

1.50 

V12 

-0.532 

4.0 

1.40 

+0.532 

4.7 

1.40 

V21 

-0.450 

4.7 

1.22 

+0.443 

5.1 

1.22 

V22 

+0.373 

3.2 

1.59 

-0.346 

3.0 

1.59 

Wll 

+0.311 

5.1 

2.25 

-0.332 

4.1 

2.43 

W12 

+0.262 

3.5 

2.62 

-0.239 

6.0 

2.71 

W21 

-0.288 

4.6 

3.09 

+0.297 

3.8 

2.53 

W22 

+0.293 

6.1 

2.43 

-0.293 

6.3 

2.62 

W31 

-0.260 

3.3 

2.25 

+0.281 

3.8 

2.34 

W32 

-0.263 

3.6 

2.62 

+0.258 

3.4 

2.43 

W41 

+0.226 

6.0 

4.40 

-0.232 

5.7 

4.21 

W42 

+0.302 

6.3 

3.28 

-0.286 

5.9 

3.37 


a Peak-peak variation in the daily mean gain, indicates range of gain drift 
during first year. 

b Mean statistical uncertainty per hour. 


Table 4. Auto-correlation Model Parameters 3 


Radiometer 

c, 

A 

B 

w(/ s pin) b 

Kll 

-0.0038 

0.0011 

0.00042 

0.966 

K12 

0.0008 

0.0011 

0.00040 

0.963 

Kail 

-0.0075 

0.0015 

0.00048 

0.960 

Kal2 

-0.0031 

0.0006 

0.00019 

0.984 

Qll 

0.0044 

0.0018 

0.00053 

0.934 

Q12 

-0.0088 

0.0007 

0.00024 

0.978 

Q21 

0.0124 

0.0088 

0.00282 

0.754 

Q22 

0.0178 

0.0128 

0.00415 

0.686 

VI 1 

0.0010 

0.0001 

0.00005 

0.989 

V12 

0.0034 

0.0014 

0.00048 

0.925 

V21 

-0.0038 

0.0010 

0.00033 

0.951 

V22 

0.0087 

0.0093 

0.00320 

0.689 

Wll 

0.0158 

0.0062 

0.00211 

0.680 

W12 

0.0048 

0.0005 

0.00019 

0.950 

W21 

0.0207 

0.0071 

0.00262 

0.644 

W22 

0.0167 

0.0053 

0.00187 

0.701 

W31 

0.0062 

0.0006 

0.00021 

0.943 

W32 

0.0077 

0.0002 

0.00007 

0.975 

W41 

0.0562 

0.0323 

0.01152 

0.374 

W42 

0.0393 

0.0194 

0.00692 

0.461 


a See equation 26 for model definition. All parameters 
are dimensionless. 

b Derived filter response at the spin frequency, 7.7 mHz. 
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a AIl values derived from a one-year simulation of WMAP data. The first 
3 data columns give the mean power in the residual map t out - from the 
simulation. The last 3 columns give an estimate of the systematic error due to 
calibration and map-making, as defined in §3.1 . For comparison, the average 


power in the CMB in each band is C 2 ~ 130 /iK 2 , (C/) 3 _ 10 


150 fx K 2 , and 



rs./ 
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Table 6. Difference Map Statistics 


DA 

Diff. map a 

C( 0 beam)/C(O) 

c 2 

QiK 2 ) 

{C;) 3 _io 

(A* K 2 ) 

( Q ) u -100 
(mK 2 ) 

|AC,| 2 b 

(MK 2 ) 

|AC;| 3 _ 10 b 

(mK 2 ) 

|AC,|i,_ioo b 

(MK 2 ) 

K1 

A12 

0.160 

107.30 

1.77 

0.12 

107.25 

1.72 

0.071 

K1 

A 34 

0.014 

6.43 

0.11 

0.06 

6.37 

0.05 

0.003 

K1 

A 1234 

0.030 

13.56 

0.30 

0.07 

13.50 

0.24 

0.011 

Kal 

Ai 2 

0.0057 

2.11 

0.14 

0.06 

2.06 

0.08 

0.002 

Kal 

A34 

0.0022 

0.11 

0.08 

0.07 

0.05 

0.01 

0.002 

Kal 

Ai 234 

0.0028 

0.01 

0.09 

0.07 

0.05 

0.03 

0.003 

Ql 

A12 

0.0035 

1.08 

0.14 

0.10 

0.98 

0.04 

0.004 

Qi 

A34 

0.0032 

0.49 

0.16 

0.13 

0.37 

0.03 

0.003 

Qi 

A 1234 

0.0044 

0.31 

0.17 

0.13 

0.19 

0.05 

0.003 

Q2 

A 12 

0.0031 

0.09 

0.14 

0.09 

0.00 

0.05 

0.005 

Q2 

A34 

0.0030 

0.28 

0.13 

0.10 

0.18 

0.03 

0.001 

Q2 

A 1234 

0.0025 

0.03 

0.11 

0.10 

0.07 

0.01 

0.004 

VI 

A12 

0.0038 

4.40 

0.35 

0.16 

4.24 

0.19 

0.006 

VI 

A34 

0.0032 

0.22 

0.35 

0.28 

0.05 

0.09 

0.016 

VI 

A 1234 

0.0024 

0.26 

0.30 

0.27 

0.01 

0.04 

0.006 

V2 

A12 

0.0043 

1.72 

0.15 

0.13 

1.59 

0.02 

0.006 

V2 

A34 

0.0026 

0.57 

0.25 

0.22 

0.35 

0.03 

0.003 

V2 

A 1234 

0.0033 

1.05 

0.38 

0.23 

0.83 

0.16 

0.008 

W1 

A12 

0.0036 

5.45 

0.36 

0.36 

5.11 

0.02 

0.025 

W1 

A34 

0.0033 

2.10 

0.74 

0.56 

1.56 

0.21 

0.019 

W1 

A 1234 

0.0034 

0.14 

0.75 

0.57 

0.39 

0.22 

0.037 

W2 

A12 

0.0029 

1.20 

0.51 

0.44 

0.79 

0.10 

0.021 

W2 

A34 

0.0026 

0.12 

0.57 

0.57 

0.44 

0.00 

0.009 

W2 

A 1234 

0.0028 

0.26 

0.56 

0.59 

0.31 

0.01 

0.025 

W3 

A12 

0.0035 

4.95 

0.49 

0.47 

4.49 

0.03 

0.013 

W3 

A34 

0.0031 

8.84 

0.90 

0.75 

8.12 

0.19 

0.039 

W3 

A 1234 

0.0039 

2.91 

1.11 

0.72 

2.20 

0.40 

0.007 

W4 

Al2 

0.0030 

1.50 

0.44 

0.48 

1.05 

0.00 

0.040 

W4 

A34 

0.0025 

0.36 

0.73 

0.64 

0.27 

0.10 

0.006 

W4 

A 1234 

0.0027 

0.73 

0.98 

0.65 

0.30 

0.35 

0.021 


difference maps from linear combinations of channels within a single DA, defined in equation (27). 
b Power in difference map in excess of white noise, |(C/) band - (G^oo-iooo I * 
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Table 7. WMAP Boresight Pointing Vectors 3 


DA/Side 

n x 

n y 

n z 

K1A 

0.0399374 

0.9244827 

-0.3791264 

KalA 

-0.0383635 

0.9254372 

-0.3769539 

Q1A 

-0.0315719 

0.9521927 

-0.3038624 

Q2A 

0.0319339 

0.9522016 

-0.3037965 

VIA 

-0.0331733 

0.9415643 

-0.3351958 

V2A 

0.0333767 

0.9414947 

-0.3353711 

W1A 

-0.0091894 

0.9394385 

-0.3425944 

W2A 

-0.0095070 

0.9458644 

-0.3244228 

W3A 

0.0098004 

0.9457678 

-0.3246956 

W4A 

0.0098081 

0.9393480 

-0.3428252 

K1B 

0.0379408 

-0.9239176 

-0.3807057 

KalB 

-0.0400217 

-0.9246344 

-0.3787473 

Q1B 

-0.0334030 

-0.9517688 

-0.3049925 

Q2B 

0.0301434 

-0.9519277 

-0.3048361 

V1B 

-0.0350363 

-0.9409454 

-0.3367405 

V2B 

0.0314445 

-0.9411385 

-0.3365553 

W1B 

-0.0114732 

-0.9388325 

-0.3441830 

W2B 

-0.0115900 

-0.9453501 

-0.3258511 

W3B 

0.0076818 

-0.9454070 

-0.3258014 

W4B 

0.0075141 

-0.9388923 

-0.3441291 


a Beam line-of-sight unit vectors in spacecraft coordi- 
nates. Available in full precision in the released time- 
ordered data. 
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Table 8. Measured Gain and Baseline Susceptibilities 3 


dg/dTwh 

Radiometer (du mK' 1 ) K _1 

db/STppA 

mKK-' 

db/dV bai 

^kv-’ 

Kll 

-0.0021 

3.52 

0.1 

K12 

-0.0185 

5.05 

3.1 

Kail 

-0.0024 

-1.47 

0.2 

Kal2 

0.0077 

2.00 

-3.2 

Qll 

-0.0037 

3.79 

-1.1 

Q12 

-0.0016 

-3.52 

-1.6 

Q21 

0.0086 

-1.00 

-2.1 

Q22 

0.0058 

-0.57 

-4.6 

Vll 

0.0018 

57.4 

32.9 

V12 

-0.0045 

-6.23 

17.2 

V21 

0.0029 

6.10 

3.4 

V22 

-0.0002 

-9.43 

-3.8 

Wll 

0.0004 

-14.7 

-43 

W12 

0.00002 

-61.9 

-11.3 

. W21 

0.0007 

-127. 

-7.4 

W22 

-0.0004 

-58.1 

0.4 

W31 

0.0003 

4.49 

5.3 

W32 

0.0021 

-20.2 

19.0 

W41 

0.0006 

41.4 

16.7 

W42 

-0.0011 

19.9 

8.0 


a The thermal values are based on fits to a 10-day cooling pe- 
riod following a partial battery cell failure. See Figure 14 and 
§3.4.1. 
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Table 9. Limits on Spin-Synchronous Environmental Effects* 


Radiometer/ 

Band 

Gain 

nK 

Thermal 

nK 

Voltage 

nK 

Kll 

-1.2 

22 

0.3 

K12 

-11.1 

32 

9.3 

Kail 

-4.2 

-9 

0.5 

Kal2 

5.8 

13 

-9.7 

QH 

-1.3 

24 

-3.3 

Q12 

-1.2 

-23 

-4.7 

Q21 

31.0 

-6 

-6.3 

Q22 

35.2 

-4 

-13.8 

Vll 

0.7 

367 

98.8 

V12 

-6.0 

-40 

51.7 

V21 

5.9 

39 

10.3 

V22 

-1.7 

-60 

-11.3 

Wll 

4.0 

-94 

-13.0 

W12 

-0.1 

-396 

-33.9 

W21 

9.1 

-812 

-22.3 

W22 

-5.8 

- 372 

1.1 

W31 

1.1 

29 

16.0 

W32 

1.1 

-129 

57.1 

W41 

21.1 

265 

50.1 

W42 

-30.1 

128 

24.0 

K 

6.2 

27 

4.8 

Ka 

0.8 

2 

4.6 

Q 

15.9 

2 

7.0 

V 

0.3 

77 

37.4 

w 

0.1 

173 

9.9 


a 1 a upper limits derived from measured gain and 
baseline susceptibilities in Table 8, combined with 
upper limits on temperature and voltage fluctua- 
tions at the spin period. Sign is preserved for each 
radiometer for roll-up by band. 
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Table 10. Upper Limits on Radiometer Cross Talk 


DA 

Electrical 

dB 

Radiometric 

dB 

K1 

-37.7 

-26.8 

Kal 

-39.5 

-30.4 

Ql 

-41.6 

-32.3 

Q2 

-41.5 

-32.2 

VI 

-43.1 

-35.2 

V2 

-42.8 

-35.4 

W1 

-48.8 

-48.3 

W2 

-47.1 

. -43.5 

W3 

-38.6 

-42.6 

W4 

-46.1 

-47.5 


I 



